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Abstract 

No terrestrial planet formation simulation completed to date has considered 
the detailed chemical composition of the planets produced. While many have 
considered possible water contents and late veneer compositions, none have ex- 
amined the bulk elemental abundances of the planets produced as an important 
check of formation models. Here we report on the first study of this type. Bulk 
elemental abundances based on disk e quilibrium studies ha ve been determined 
for the simulated terrestrial planets of O'Brien et al.l ( 20061 ). These abundances 



are in ex cellent agreement wit h observed planetary values, indicating that the 
models of O'Brien et al.l ( 2006f ) are successfully producing planets comparable to 



those of the Solar System in terms of both their dynamical and chemical prop- 
erties. Significant amounts of water are accreted in the present simulations, 
implying that the terrestrial planets form "wet" and do not need significant wa- 
ter delivery from other sources. Under the assumption of equilibrium controlled 
chemistry, the biogenic species N and C still need to be delivered to the Earth 
as they are not accreted in significant proportions during the formation process. 
Negligible solar photospheric pollution is produced by the planetary formation 
process. Assuming similar levels of pollution in other planetary systems, this 
in turn implies that the high metallicity trend observed in extrasolar planetary 
systems is in fact primordial. 
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1. Introduction 

Terrestrial planet formation, both in terms of the dynamics and chemistry 
involved, is still not fully understood. Dynamic ally, basic p l aneta ry formation is 
described through the planetesimal theory (see Chamberg ( 2004 ) for a detailed 



review). This theory sees planetary formation occurring through three main 
steps. Initially, dust settles into the midplane and accretes to form planetesi- 
mals, the first solid bodies of the system. This stage is followed by the coUi- 
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sional accretion of planetesimals to produce planetary embryos. The growth of 
the embryos occurs initially via runaway growth (where an increasing geometric 
cross-section and gravitational field allow for the accretion of an ever increas- 
ing number of planetesimals) before transitioning to oligarchic growth (where 
neighboring embryos grow at similar rates). Finally, as numbers of planetesi- 
mals decrease, the interaction between embryos becomes the dominant factor 
as they perturb each other onto crossing orbits, thus producing accretion via 
violent collisions. 

Many attempts have been made to simulate the third stage of terrestrial 

plane t formation descr ibed above (e.g. iKominami and Idal200^ 2004 ; Chambers and Wetherill 
19981: IChambCT^ I2OOII : Ravmond et al.F2004l ). However, these simulations have 



had limited success and no simulation has been able to exactly reproduce the ter- 
restrial planets of the Solar System in terms of their number, masses and orbital 
parameters. For example, direct N-body simulations have produced systems 
with approximately the correct number of planets y et with an orbital excita- 
tion greater than that observed for the Solar S ystem (jChambers and Wetherilll . 
1998t IChamberd . l200lt [Raymond erall. |2004) Similarlv. simulations mcorpo- 



rating tidal torques (e.g. Kominami and Ida l2002i [2004 ) have produced too 
many planets but with more favorable excitation levels. 

Recent developments in modeling have occurred with the incorporation of 
dynamical friction, the process whereby equipartitioning of energy between low 
mass planetesimals and larger embryos results in reduced relative velocities for 
the embryos, thus increasing their probability of accreting. Dynamical friction 
has been shown to be a viable mechanism to reduce the dynamical excitation 
levels of the final planets to bet ter agree with observed values within the So- 
lar System ( Levison et al. . 20051 ). Several authors have recently performed N- 
body simulations with a high enough resolution to accurately m odel the effects 
of dyn ami cal friction during terrestrial plane t accretion, such as 
' 2006|) and lRavmond et all (|2006l . [2007ll2009t ). The simulations of 



O'Brien et al 



O'Brien et al 



20061 ) ■ which will be used in this work, represent a factor of ~5 increase in the 



number of gravitationally in teracting bodies compared to most other previous 
simulations of this type (e.g. Chamber^ 200 1 ) . and because of the large number 
of small planetesimals in the simulation, dynamical friction is more accurately 
treate d than in thos e previous simulations. The terrestrial planets produced 
in the lO'Brien et al.l simulations are a significantly better fit to the observed 
properties of the terrestrial planets of the Solar System than those of previous 
simulations. There are fewer planets (averaging 3.25 planets per simulation) 
and the dynamical excitation of the planetary systems is comparable to that of 
the actual terrestrial planets in the Solar System. Furthermore, the accretion 
timescales of the planets simulated are in agreement with the acc retion timescale 
for th e Earth (~60 Myr) as obtained from i^^Hf-^^^W dat ing (jTouboul et al 



2003). The simulations o flRavmond et all (|2006ll2007il2009l ). which have a com- 
parable resolution to the lO'Brien et al.l simulations, also produce a comparable 
improvement over previous simulations. Thus it can be seen that these recent 
N-body simulations have provided us with some of the most realistic and feasible 
models of terrestrial planet formation completed to date, and clearly represent a 
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significant step towards understanding the last stage of the planetary formation. 

Several other types of terrestrial planet formation models have been devel- 
oped recently as well. The effects of the depleting Solar nebula gas on terres- 
trial planet accretion, both in terms of its effects through gravitational tides 
and also by causing the sweeping of secular res onances through t he te rrestrial 
planet r egion, has been s t udied analytically by iNagasawa et aL ( 2005 ). More 
recentlv. iThommes et al. ( 2008 ) incorporated these effects into an N-body in- 
tegrator. That work did not incorporate a planetesimal swarm (and hence did 
not include dynamical friction), but found that under certain circumstances, the 
masses and dyn amical state of the terrestrial plan ets in the Solar System could 
be reproduced. iBromlev and Kenvon ( 20061 ) and Kenvon and Bromley (l2006h 
developed a hybrid code in which large bodies are treated via an N-body inte- 
gration scheme, while small bodies are treated with a multi-annulus coagulation 
model that can treat the growth of bodies starting from km-scale planetesimals. 
They too found that their model could reproduce many of the dynamical char- 
acteristics of the terrestrial planets in the Solar System, and while the final 
eccentricities of the planets were somewhat high, it is likely that a more com- 
plete treatment of the effects of fragmentation, and of tidal drag on large bodies 
from the nebular gas, would help to decrease them. A full compari son of these 
models with the direct N-body simulations of lO'Brien et al.l (|200d ) and others 
is beyond the scope of this article, but we will discuss the implications of those 
models for the work presented here in Sec. H) 

While models have advanced significantly in their ability to reproduce the 
terrestrial planets from a dynamical point of view, the vast majority of terrestrial 
planet formation studies completed to date have been limited (at best) in their 
simultaneous chemical composition studies. Several studies have estimated the 
amount of potential ly hydrated material that may have accreted into the ter- 



restri al plan ets (e.g. Rayniond et al. 2004 O'Brien et al. 2006t Raymond et al 



2007[) while lO'Brienet all (|2006f) and [Raymond et al.l \20m also considered 



the delivery of siderophile-element-rich late veneer material need ed to accoun t 
for the highly siderophile element budget of the Earth's mantle (jChoul . 119781 ). 
However, detailed examination of the bulk chemical composition of the resulting 
planets as a test of dynamical simulations has never been thoroughly explored. 
This combination of the two approaches to produce a comprehensive model of 
both the chemistry and dynamics of terrestrial planet formation is essential to 
determine how well current numerical simulations reproduce not only the masses 
and dynamic state of the terrestrial planets, but also their bulk composition. 
The present study represents a first step towards such a model. 

We have derived detailed bulk e lemental abundances for all of the terrestrial 
planets formed in the simulations o Brien et al. I (|2006l) . These predicted com- 
positions are compared to the bulk compositions of the actual terrestrial planets 
as an examination of the chemical plausibility of the dynamical simulations cur- 
rently being used to model terrestrial planet formation. This approach allows 
us to not only examine the bulk elemental composition of the final planets but 
to also study the compositional evolution with time. Additionally, we also in- 
vestigate the delivery of hydrated material, which has important implications 
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for the development of habitable terrestrial planets, along with the composition 
of the "late veneer" material accreted by the planets. Finally, the amount of 
material accreted by the Sun as "pollution" during the terrestrial planet for- 
mation process and the resulting changes in stellar photospheric abundances 
are also examined. This study is the first to attempt to produce a compre- 
hensive, simultaneous model of both the chemistry and dynamics of terrestrial 
planet formation. As such, it represents a significant step forward in producing 
a cohesive model of terrestrial planet formation. 



2. Simulations 



2.1. Dynamical 



The eight high-resolution N-body accretion simulations of lO'Brien et al 



(2006) were utilize d in this study. These simulations were performed with the 
SyMBA integrator ( Duncan et al.l . [l998h . which is ideally suited for this type of 
simulation as it is designed to rapidly and accurately integrate the orbits of a 
system of bodies on Keplerian orbits, and can treat close encounters and merg- 
ers between colliding bodies. Each simulation consists of an equal distribution 
of mass between Mars- mass embryos (O.OOSSM^, 25 per simulation) and plan- 

etesimals l/40th the size of the embryos (O.OO233M0, ~1000 per simulation). 
As is common in this type of simulation, the planetesimals gravitationally inter- 
act with the embryos, but not one-another. The initial distribution of bodies is 
constrained by the disk surfa ce density profile X^l'")— So(^'/lAU)^'^/^ with J2o 
= 8 gcm~^ (Chamber^. 2001 ) and bodies initially located between 0.3 and 4.0 
AU ( O'Brien et al.l 20061 ). Four simulations were run with Jupiter and Saturn in 
low-eccentricity, low-inclination orbits slight ly closer together than their present 
locat i ons, as predicted by th e Nice Model (jGomes et al. . l2005t iTsiganis et al 



boost iMorbidelli et al.l . l2005l ) (hereafter termed 'Circular Jupiter and Saturn', 
CJS) and four with Jupiter and Saturn in their current, slightly eccentric or- 
bits (hereafter termed 'Eccentric Jupiter and Saturn', EJS). Multiple terrestrial 
planets formed within 250Myr in all simulations. The general system structure 
is shown in Figure [1] and the planetary orbital properties are shown in Table [TJ 
It is currently unclear whether CJS or the EJS (or some other c onfiguration 



best r e presents the in i tial c onfigu ration of Jupiter and S aturn. I Gomes et al 
(|2005h . lTsiganis et all (|2005l ). and iMorbidelh et all (|2005h have shown how the 
CJS configuration of the giant planets, and their subsequent orbital evolution 
via interactions with the primordial trans-Neptunian disk, can explain their cur- 
rent orbital configuration, the populations of Jupiter T rojan asteroids, and the 
Late Heavy Bombardment of the inner Solar System. IMorbidelli et al. ( 2007 ) 
has shown that such a configuration is a natural consequence of the giant plan- 
ets' formation in, and interactions with, the protoplanetary disk. However, 
simulations of terrestrial planet formation tend to find that some aspects of the 
terrestrial planets, namely a compact system and a small Mars, are best repro- 
duced with a configuration more like in the EJS simulations, wh ere Jupiter has 
an eccentricity comparable to or larger than its current value (|0'Brien et al.L 
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2006t iThommes et al.l . bOOS i: iRavmond et al.l . [200i). Given the current uncer- 
tainty, we study both the CJS and EJS simulations here, as they give a good 
representation of the range of possible initial configurations for Jupiter and 
Saturn. 

It should be noted that t he formation o f Merc ury analogs is not currently 
possible in the simulations of lO'Brien et al. I (I2OO6I) as all of the embryos begin 
with a mass almost twice that of Mercury. Additionally, the planetesimals 
are not gravitationally interacting with each other, preventing formation of a 
small 'Mercury' via planetesimal accretion. Finally, perfect accretion is assumed 
throughout the simulations, preventing giant impact events from striping off 
the crust and mantle of a differentiated proto-Mercury, a scenario that widely 
suggested to explain Mercury's current high density. As the simulations, results 
and their implications have already been discussed in gr eat detail by previous 
publications, the reader is referred to lO'Brien et al.l (2000) for further discussion 
of these aspects. 



2.2. Chemical 

The chemical composition of material within the disk is assumed to be de- 
termined by equilibrium condensation within the primordial solar nebula. This 
assumption is justified as numerous detailed analyses of primitive chondritic 
material have clearly shown that their bulk elemental abundances are smooth 
functions of their equilibrium conde nsation temp erature, as determined for a so- 
lar composition gas at low pressure (iDavid . I2OO6I ). Several possible explanations 
for this pattern have been proposed, with the most widely accepted being that 
the chemistry of their paren t planetesimal s was established by midplane temper- 
ature and pressure profiles ( CassenL 2001). Alternatively, the sluggish nature of 
reaction kinetics within the cooler regions of the outer nebula may also produce 
the observed trend. Although the exact cause of the primitive chondritic pattern 
is still unclear, we are able to utilize the fact that equilibrium condensation tem- 
peratures are reliable proxies for the bulk elemental abundances of rock-forming 
elements within the early Solar System. Furthermore, observational evidence 
of these equilibrium compositions is still seen today, preserv ed as the thermal 
stratification of the asteroid belt ( Gradie and Tedesc ol ll982h . The mineralogy 
of primitive chondritic meteo rites is simi lar in most respects to that predicted 
by equilibrium condensation ( Ebell . 2006[ ). Thus the assumption of equilibrium 
composition is a valid starting point for this type of study. This assumption 
in turn implies that the primary controls on embryo and planetesimal composi- 
tions are the radial midplane temperature and pressure gradients encountered 
by the nebular material. 

In order to determine the equilibrium composition of the solid material, we 
utilized the commercial software package HSC Chemistry (v. 5.1). HSC Chem- 
istry determines the equilibrium chemical composition of a system by iteratively 
minimizing the system's Gibbs free energy, using the GIBBS equilibrium solver, 
as described in White et al.l ( 19581 ). Thi s software has been successfully used in 
recent studies of solar nebula chemistry (jPasek et al.l . Eooi) and supernova stel- 
lar outflows, producing compositions that correlate with observed mineralogy in 
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presolar grains found in interplanetary dust particles (jMessenger et all 120051 ). 



As we are working with similar spatial, thermal and pressure scales as these pre- 
vious studies and as the key theory applied here is not spatially dependent, we 
feel confident in applying it to this current study. Several limitations have been 
encountered with the program, primarily the limited number of elements and 
species possible to incorporate in a simulation and its inability to handle species 
in non-ideal solid solutions such as pyroxene and olivine, two of the major rock- 
forming minerals on the terrestrial planets. End member species can and have 
been considered but the non-ideal solid solution interaction between them is not 
currently possible. These limitations will not significantly affect the final con- 
clusions of this current study but will be the subject of future work. The list of 
solid and gaseous species included in the HSC Chemistry calculations are shown 
online in Table [7] (online supplemental material). Note t hat although silicate 



hquid s were found to be stable in dust-enriched systems (jEbel and Grossman , 
200Cl[ ) , here we are only considering compositions with a solar dust/gas ratio E ■ 



As such, liquids do not play an important role in the condensation sequence and 
are thus neglected from the current calculations. 

Our present simulations incorporate 14 major rock-forming elements (C, N, 
O, Na, Mg, Al, Si, P, S, Ca, Ti, Cr, Fe and Ni), along with H and He, which dom- 
inate the partial pressures of protoplanetary disks. Current solar photospheric 
abundances are utilized as a proxy for the initial composition of the solar nebula. 
All inputs are assumed to initially be in their gaseous, elemental forms. Sol ar el- 
emental abundances for all elements were taken from Asplund et al. ( 20051 ) and 
are shown in Table [21 Although there h as been some recent discussion a bout 
revising the Solar C and O abundances ( Socas-Navarro and Norton . 20071 ). we 
have adopted the most recent and widely accepted values for this current study. 
The 50% condensation temperature calculated for each element in a gas with 
a Solar photosphere composition is shown in Table [S] The 50% condensation 
temperature is the temperature at which 50% of the original element has con- 
densed, with the remaining 50% still in the gaseous phase. As such, it provides 
an effective measure of the volatility of a given element in the system. It should 
be noted that although O is present in high temperature condensates and sil- 
icate species, it does not obtain 50% condensation until water ice condenses 
within the system, hence its low 50% condensation te mperature . Add itionally, 
the 50% condensation temperatures from the models of Lodders J 20031) are also 
provided. The average difference between this study and that of Lodderd ( 20031 ) 
is just -1-11. 75K (defined as Tprcsent study - TLoddcrs), with the largest difference 
being 83K for P. The e xcellent agreem ent between the temperatures obtained in 
this study and those of lLodder^ ( 20031 ) implies that the chemical models utilized 
in this study are accurately reproducing the thcrmodynamical condensation se- 



^Here the dust/gas ratio is taken to be the cosmochemical definition, that is the atomic 
Si/H ratio of 3.5xl0~^. The usual astrophysical definition of the dust/gas ratio is the ratio 
of mass of the solid material to the mass of the gas within the system. As expected, this value 
varies over the disk as condensation progresses. 
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quence of the solar disk. 

In order to relate the chemical abundances to a spatial location within 
the original disk, we used midplane pressure and temperature va l ues o btained 



from the axisymetr ic a viscosity disk model of iHersant et al.l (|200lh . The 



Hersant et alj (|200lD model is a two-dimensional time-dependent turbulent ac- 
cretion disk model incorporating vertical disk structure, turbulent pressure and 
self-gravity. As for other standard disk models, a Keplerian rotation law, hydro- 
static equilibrium and energy b alance between visco us heating and cooling due 



to radiative losses are assumed ( Hersant et al.l . 120011) . The effects of irradiation 



from the central star are neglected, as are other disk features such as inner holes 
and shadow zones. In the current study, we limit the midplane conditions to 
be those produced by the "nominal" model defined bv lHersant et al. ( 2001 ) to 



be a stellar mass accretion rate of 5x10 ^M0yr ^, an initial disk radius of 17 

AU and an a value of 0.009. Previous work has successfully applied this model 
to constrain the bulk composition of Jupiter (e.g. Pasek et al.i]2005i ). Pressure 
and temperature values were determined with an average radial separation of 
0.03AU throughout the study region. Due to computational limitations, pres- 
sure values were rounded to the nearest order of magnitude. This is not expected 
to greatly alter the final planetary compositions produced as temperature exerts 
a far greater control on the equilibrium composition at a given radii. 

The midplane temperature and the pressure, and thus also the equilibrium 
composition of solids present within the disk, changes with time as the disk 
evolves. In order to capture this effect in our simulations, we determined an 
ensemble of predicted planetary compositions constrained by disk conditions 
at multiple time periods. At the earliest, we use the temperature and pressure 
profiles obtained for disk conditions at t = 2.5 x lO^yr, where conditions were first 
determined to be suitable for solids to be present across the entire radial region 
being modelled in the dynamic simulations. We end with disk c onditions at t 
= 3x lO^yr, the average lifetime for the protoplanetary gas disks jHaisch et al 



20011) . Five cases between these end points are considered (for disk conditions at 
t = 5xlOV, IxlO^yr, 1.5xlOV, 2xl0^yr, 2.5xlOV)- Thus we determine the 
planetesimal and embryo compositions for a total of seven cases, covering the 
entire range of times during which embryos and planetesimals could potentially 
form, either early or late in the lifetime of the disk. Note that the timescale used 
for the evolution of the disk conditions is not coupled with the timescale of the 
dynamical simulations. Rather, chemical compositions are simply determined 
for a range of disk mid-plane conditions and no time- variation in disk conditions 
over the duration of the dynamical simulation is incorporated. 

2.3. Combining Dynamics and Chemistry 

In order to combine these two different modeling approaches, we assume 
that each planetary embryo and planetesimal of the dynamical models retains 
the chemical composition in equilibrium with the nebula in the region that it 
first formed. The bulk compositions of the final planets are simply the sum of 
each object they accrete. By tracing the origi n of each emb r yo an d planetesi- 
mal incorporated into the final planets of the O'Brien et al.l ( 20061 ) dynamical 
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simulations, and calculating the chemical composition of those bodies based on 
their original locations, we constrain the bulk composition of the final terres- 
trial planets. This procedure was completed for all planets at each of the seven 
different time steps simulated within the chemical models. 

2.4- Stellar Pollution 

Stellar pollution refers to the late addition of material to the stellar photo- 
sphere, usually during planet formation and/or migration. It has been suggested 
as being a possible cause for the high mctall icity trend observed in known extra- 



solar planetary host stars ( Laughlinl . I2OO0I : iGonzalez et al.l . l200lHMurrav et al 



200lh . The ex tent of stellar poUutfo n produced by terrestrial planet formation 
simulations of lO'Brien et al. ( 20061 ) was measured by determining the amount 
and composition of material accreted by the Sun during the formation process 
and the resulting photospheric elemental abundance changes such an addition 
would produce. Any solid material migrating to within 0. 1 AU from the Sun is as- 
sumed to have accreted onto the Solar photosphere through gravitational attrac- 
tion. This material is then assumed to have been uniformly mixed throughout 
the solar photosphere and convective zone. Gravitational settling and turbu- 
lence within the convective zone are neglected from this present study, primarily 
due to the large time scales involved in settling (e.g. ^10% reduction in photo- 
spher ic heavy metal abundances by the current solar age ( Pinsonneault et all . 

and the uncertainties involved in the effects of turbulence on specific 
elements. Such models are beyond the scope of this paper. However, it is 
expected that each of these processes would reduce the amount of accreted ma- 
terial present within the upper layers of the Sun, thus reducing the photospheric 
elemental abundance. As such, the resulting abundance changes determined in 
this study represent the maximum possible values that we would expect to ob- 
serve in the Solar spectrum. Additionally, at earlier times, the mass of the solar 
convective zone would have been greater than the current value of 0.03 
again resulting in a smaller change in photospheric abundance than predicted 
by this current study. 

The mass of each element accreted by the Sun was determined in the same 
way as described in section [2?3l for terrestrial planets. The resulting photospheric 
abundance is given by: 



[X/H] = log 



fx 



(1) 



where [X/H] is the resulting abundance of element X after accretion of ter- 
restrial planet material, fx is the mass abundance of element X in the Solar 
photosphere after accretion and f^^ q is the initial mass abundance of element 

X in the Sun before accretion ffrom [Murray et "aL 2001 ). fx (7) values were 

obtained by utilizing the solar abundances o flAsplund et al. ( 2005 ) and a cur- 
rent solar convective zone mass of 0.03Mq (iMurrav et all . l200l[ ). The present 
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approach only addresses pollution by the direct accretion of planetesimals and 
embryos during the current simulations, yet pollution may also occur both be- 
fore planetesimal accretion (and thus the simulation) begins and through the 
accretion of dust produced by the impact and coUisional events during the for- 
mation process. Pollution by these processes is believed to be small, resulting 
in at most a factor of two increase in the amount of material added to the solar 
photosphere. 



3. Results 

3.1. Abundance Trends 

The bulk elemental abundances for all simulated terrestrial planets (as ele- 
mental wt% of the planet) for each set of disk conditions examined are provided 
online as supplemental material (Table [8|). Additionally, the Mg/Si value for 
each planet is also listed. 

From these results, it can be seen that the bulk elemental abundances of 
the final terrestrial planets vary significantly with disk conditions, especially for 
certain elements such as O, Al and Ca. As expected, the relative abundances 
of the most refractory elements (Al, Ca, Ti) are observed to decrease for disk 
conditions at later times while the more volatile elements (O, Na, H) all increase 
in abundance. The largest variation occurs in the O abundance, with simulations 
for disk conditions at t = 3x lO^yr producing terrestrial planets with up to 30.79 
wt% more O than identical simulations for disk conditions at t = 2.5xl0^yr. 

These compositional variations can be directly attributed to the changing 
conditions of the disk itself. As the temperature and pressure at the mid-plane 
decrease with time, the equilibrium composition at a specific radius also changes. 
Thus the assumed composition of the planetesimals and embryos also change, 
in turn producing variations in the final planetary abundances. Consequently, 
as time progresses and the disk cools, the planetary composition can be seen to 
decrease in abundance of the most refractory elements (Al, Ca) and increase in 
the more volatile elements (O, H) as volatile, hydrated species become stable 
throughout an increasingly larger fraction of the disk. 

A more accurate picture of the predicted planetary abundances is obtained 
by examining the abundances produced for each specific set of disk conditions 
simulated. Figure [2] shows the Si and planet-normalized bulk elemental com- 
position for the planets produced in the CJS-1 and EJS-1 simulations for all 
seven disk conditions. The simulated planets were normalized to the bulk ele- 
mental composition of Venus, Earth or Mars on the basis of their orbital prop- 
erties, primarily their se mimajor axis. Reference So lar Syste m planetary abun- 
dance s were taken from Morgan and Anders ( 19801 ) fVenus). iKargel and Lewis! 



("igga") (Earth) and lLodders and FeglevI (|l997t ) (Mars). The normalized values 
were obtained via: 
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/ Wt. %Qf X '\ 

V wt. % of Si 



TVT Till r 1 , ^ ^ simulated planet 

Normalized abundance or element A ~ r (2) 

'wt.%ofX^ 



wt. %of Si 



terrestrial planet 



Values above unity indicate an enrichment in the simulated terrestrial planet, 
while those below it indicate depletion. 

All three suites of planetary reference values are themselves based (to some 
degree) on simulations. Bulk Earth values were obtained by extrapolation of 
Best Bulk Silicate Earth (Best BSE or BBSE) abundance values (as obtained 
from previously published studies of mantle xenoliths, pyrolite and basaltic el- 
emental ratios), the observed volatility trend for volatile lithophile elements 
and the known physical and seismic properties of the interior of the Earth 
( Kargel and Lewisl . ll993l) . This approach, although necessary, does introduce 
uncertainties especially when calculating the composition of the core as par- 
tioning ratios are not well known for all elements. Furthermore, the use of a 
volatility trend requires knowledge of the condensation temperature of the el- 
ements. While such information is well known for several species (such as Fe, 
Mg and Si with errors on the order of 3% of the condensation temperature), 
larger degrees of uncertainty remain for more volatile elements. This uncer- 
tainty produces errors in the predicted condensation temperatures of ±14% of 
the condensation temperature itself and results in final elemental abundance er- 
rors of 5-10% of the actual abundance for most elements, reducing to less than 
1% for Fe, Mg and Si. Of the elements considered in the current study, those 
with the greatest overall uncertainty are O, S and C. Although not quantified by 
Kargel and Lewii ( 1993 ). the uncertainty in O is derived from the fact that the 



O abund ance was obtained based on the valence states of cations. The approach 
taken bv lKargel and LewisI ([l99i) in obtaining their bulk compositional values 
is predicated upon the assumption that 15% of all Fe present is in the form 
Fe'^^ . As such, any variation within the assumed valence state distribution will 
induce a variation in the O abundance of Earth. C and S abundances similarly 
suffer from large uncertainties (~25% of the given abundance) due to variations 
in the published BSE values, the predicted condensation temperatures and the 
core partitioning coefficients. 

Bulk Mars abundances were taken from iLodders and Fegleyl (|l997l ) and are 
based on the combination of known meteorite material (H, CV and CI meteoritic 
material) required to reproduce the oxygen isotopic abundances observed for 
the SNC meteorites. Variations in the elemental compositions within each class 
of meteorite are thus the main source of uncertainty in the current estimates. 



Lodders and FeglevI (|1997( ) estimate their final abundance errors to be ±10% 



for all elements studied. 

Finally, the marked lack of data regarding the elemental abundances of 
Venus resulted in the necessary use of t he purely theore tical bulk planetary 
abundances of iMorgan and Anders! (|l98nf ) for this study. I Morgan and Anders! 
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(1983) adopt a similar approach to the one used here and assume that the com- 
position of soHd material initially present within the system is controlled by 
equilibrium condensation. The exact predict ed composition was o btained from 
the Ganapathy- Anders 7-component model ( Morgan et al. . 19781) . This model 



determines the composition of a body based on the amount of early condensate 
(i.e. first solids present within the disk), metal, silicate, troilite, FeO, MnO and 
Cr203 believed to be incorporated into a final planetary body. Bulk elemental 
abundances are thus inferred from general cosmic proportions as compared to 
an "index element" for each component. As such, general and solar elemen- 
tal ratios are assumed to be homogenou s and remain consta n t for the entire 
planetary system. For the simulations of Morgan and Anderd (1983), the lim- 
ited information available on the index elements for Venus was supplemented 
by geochemical restrictions from the Earth, Moon and chondritic meteorites 
to obtain predicted bulk planet abundances. As such, large uncertainties exist 
in all elements studi e d. As quantified uncertainty factors are not provided by 
Morgan and Anderd (|l98d) . it is impossible to gauge the errors on their pre- 



dicted abundances. However, they are believed to be considerable, resulting in 
elemental abundances for Venus that should be taken as a guide only. 

From Figure[2]it can be seen that although no single simulation exactly repro- 
duced the compositions of the terrestrial planets, the time when the adopted disk 
conditions produced simulated compositions that most closely resemble those 
of the terrestrial planets (the 'best fit' time) is SxlO^yr, based on the compar- 
atively small deviations produced in the bulk elemental abundances. For disk 
conditions at later times, significantly larger enrichments can be seen in O, Na 
and S for all planets. Although not shown here, disk conditions at this same time 
also produces the best agreement between simulated and observed abundances 
for the other six simulations studied and as such future discussions will focus 
on the compositions produced by the disk conditions at this time. This result 
does not imply that the material which formed the terrestrial planets actually 
condensed out of the Solar nebula exclusively at 5xl0^yr. Rather it is simply 
the time at which the disk conditions and resultant snapshot of the chemistry of 
the Solar disk (as used here) best reproduced the expected abundances. Similar 
plots for the other six simulations studied are included as supplemental material 
online. 

The elemental wt% produced for each of the simulated planets based upon 
disk conditions at 5x lO^yr are shown in Tabled) Numerically, these abundances 
agree with Solar System values to within 1 wt% for all elements studied except 
for Mg (up to 2 wt% deviation), Fe and O (up to 2.5 wt% deviation) and S (up to 
5 wt% deviation). In terms of normalized values, the predicted abundances are 
within a factor of 0.5 for almost all elements studied. Enrichments can be seen 
in Na (up to 4x the expected value) and S (up to 8x the expected value) and 
are discussed below. These abundances match well with observed Solar System 
values, indicating that these terrestrial planet formation models are producing 
planets with both orbital properties and bulk elemental abundances comparable 
to those of the terrestrial planets of the Solar System. 

Despite the excellent agreement between the observed and predicted abun- 
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dances, the planet-normalized spider plots of Figure [H clearly display several 
systematic deviations from the expected terrestrial planet abundances. Specif- 
ically, considerable deviations can be seen in P (Venus analogs only), Na and 
S for all simulations. This is partly due to the uncertainties in the expected 
planetary abundances for these elements as was previously discussed. Addition- 
ally, as P, Na and S are some of the most volatile elements in this system, it 
is believed that the observed enrichments are a result of the fact that we do 
not consider volatile loss during the accretion process. Both the dynamical and 
chemical simulations assume that in any given impact, perfect merging of the 
two bodies occurs and all mass is retained. However, such violent impacts are 
known to both melt and eject a considerable portion of both the target bod; 



and the impactor (such as in the moon forming impact (jHartmann and Davii 



I 



19751: ICameron and Wardl . \l97^ ). This has the potential to drastically alter 



the bulk elemental abundances of a planet. It is expected that a significant 
fraction of the most volatile species in the system will be lost from the planet 
to the surrounding disk during the accretionary process, thus reducing the ob- 
served enrichments. This effect has not been incorporated into our simulations, 
resulting in volatile enriched final planetary bodies. The results of a first or- 
der approach to include the abundance effects of volatile loss are discussed in 
Section [S^Hl 

Some degree of radial compositional variation is captured by the current 
models. Simulated planets deemed to be Venus and Earth analogs were found to 
produce equivalent normalized abundances for the major rock forming elements 
when normalized to the elemental abundances of Venus and Earth. Differences 
in the quality of the fit occur for P, Na and S and are due to uncertainties 
in the abundances of Venus. The Mars analogs are found to have a slightly 
better fit when normalized to the elemental abundances of Mars, as compared 
to those of Venus and Earth. This difference is minor with Earth abundances 
arguably producing an equivalent normalization. However, given that the radial 
compositional differences observed for Venus, Earth and Mars are small and that 
planetary accretion is stochastic and includes a significant amount of radial 
mixing of material, it is not surprising that the approach adopted in this study 
has been unable to complectly reproduce the expected radial variations. 

3.2. Planetary Geochemical Ratios 

In addition to matching the bulk elemental abundances, simulations such as 
these should also reproduce the key geochemical ratios for the terrestrial planets. 
Here we have considered the values of Mg/Si, Al/Si and Ca/Si. All three of these 
values differ by chemically significant amounts from the expected planetary 
values for all planets produced. For simulations run with disk conditions at 
times of 5x10^ years and beyond, the planetary ratios are identical to the solar 
input ratios. This effect is produced by the fact that for disk conditions at t = 
5x10^ years, all the Mg, Al, Ca and Si has condensed out over the study region. 
This results in the solid species (and thus also the planets produced) possessing 
the same chemical ratio s as present in the initial solar nebula. For example, 



for lHersant et al.l (|200lh disk conditions at t = 5x10^ years, Mg/Si, Al/Si and 
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Ca/Si ratios in the solid material reach Solar values by 0.5AU, clearly interior to 
the primary feeding zone for the formation of Earth and Mars. Unfortunately, 
both the Mg/Si and Al/Si values are lower than is observed for the Earth but 
above current Martian values while the Ca/Si values are above those of both 
Earth and Mars (but in agreement for Venus) (see Figure [3]). Thus for the 
majority of disk conditions presently considered we are producing planets that 
have compositional ratios in between those of Earth and Mars and are slightly 
enriched in Ca. 

Variations in these elemental ratios from solar ratios are produced in simula- 
tions based on disk conditions at 2.5x10^ years (see Figure[3]). At this time, the 
inner disk (within approximately lAU) is dominated by Al and Ca rich species 
such as spinel (MgAl204) and gehlenite (Ca2Al2Si07), resulting in planets with 
high Al/Si and Ca/Si values and low Mg/Si values as can be seen in Figure [3l 
Thus it can be seen that there is a temporal variation in the chemical ratios of 
the planets produced. This further supports our choice of disk conditions at t 
= 5x10^ years as producing the "best fit" abundances as although the ratios 
discussed here do not precisely agree with the observed values, the deviation is 
significantly smaller than for earlier disk conditions. 

The fact that the simulated chemical ratios do not agree with observed plan- 
etary values suggests two possible causes. The first is that the dynamical sim- 
ulations are not forming planets from material sourced from the correct region 
of the disk. There is a small radial region within the disk where the values 
of Mg/Si, Al/Si and Ca/Si are in agreement with those of Earth. This region 
occurs between 0.61 to 0.68AU (corresponding to temperatures of 1352 and 
1305K) for disk conditions at t = 2.5x10^ years a,nd 0.12 to 0.13AU fo r disk 
conditions at t = 3x10^ years for the disk model of Hersant et al.l ( 200lh . Thus 
one possible solution for the current difference may be that the dynamical mod- 
els need to form Earth from material located within this region. However, it is 
exceedingly difficult to imagine how such a formation process would occur as it 
requires the movement of a large amount of material over a significant distance 
through the disk. Current formation simulations do not produce this degree of 
radial mixing and furthermore it is questionable whether sufficient material to 
form the terrestrial planets would even be located within this region, making it 
unlikely that such a scenario is feasible. Additionally, bulk elemental planetary 
compositions produced at these temperatures would not be in agreement with 
observed bulk planetary elemental abundances. Finally, observations have sug- 
gested that t he temperature at l AU in the disk of young stellar objects is less 
than 400K (.Beckwith et al.l . ll99(J) . far less than the 1300K required to produce 
the necessary Mg/Si, Al/Si and Ca/Si values at lAU. This suggests that the 
dynamical models utilized here are not the source of the deviation in Mg/Si, 
Al/Si and Ca/Si values. 

The second possibility is that the composition of solids within the disk are 
controlled by disequilibrium condensation. Disequilibrium condensation refers 
to the chemical model in which once a solid has condensed, it is removed from the 
system and may no longer interact with the remaining gas. Such a process may 
occur if condensation occurs rapidly enough to for large bodies to grow quickly 
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and thus shield the interiors from further equihbrium reactions (jCowlevl . Il995f ). 
Under s uch coridition s, secondary condensates would have Mg/Si values above 
average ( Searsl 2004 ). thus possibly increasing th e Mg/Si value t o the required 



Earth value over the primary feeding zone in the O'Brien et alj ( 20061 ) simula- 



tions. The possible role of disequilibrium processes does not invalidate our initial 
assumption of equilibrium controlled abundances. The assumption of equilib- 
rium is acceptable for determining bulk elemental abundances but finer details 
of planetary composition will clearly need to incorporate a greater interaction 
between a variety of equilibriums and disequilibrium processes. Furthermore, 
equilibrium-based predictions for simulations of this type need to be calculated 
initially to provide a baseline for future disequilibrium studies. 

3.3. Oxidation State 

The oxidation state of the planets is of great chemical and geological in- 
terest. As we are currently calculating predicted bulk elemental abundances 
and not detailed mineralogies for the simulated planets, estimates of the bulk 
oxidation state of the planet are obtained by calculating the oxidation state of 
the embryos and planetesimals before accretion occurred based on their equilib- 
rium compositions. The resultant oxidation states are shown in Figure H) The 
current simulations are producing very reduced planetary compositions for the 
earliest set of disk conditions studied, containing large amounts of metallic Fe 
and little FeO. For disk conditions at later times, the compositions become more 
oxidized, evolving through the oxidation state of the H-type meteorites to finish 
closer to the redox state of the CR meteorites, primarily due to the increased 
amount of magnetite (Fe304) and fayalite (Fe2Si04) accreted by the planets. 
For midplane conditions at t = 2.5x10^ years, fayalite is only present in signif- 
icant amounts beyond 1.07AU while magnetite is only present beyond 4.82AU, 
beyond the dynamical simulation range. For conditions at t = 3x10^ years, 
fayalite is present beyond 0.2AU with magnetite now present beyond 0.95AU. 
Thus for later disk conditions, the feeding zones of the terrestrial planets (and 
thus the planets themselves) are more oxidized. 

In addition to the observed temporal variations seen in the current model, 
there is some evidence to support the idea t hat the localized oxida tion state var- 



ied throughout the solar disk. For example. ISchrader et al.l (|2008f ) outlines how 
the type I chondrules present in CR chondrites formed in a reducing environment 
while the type II chondrules (also present in CR chondrites) formed later under 
oxidizing conditions. However, these variations in the ambient oxidation state 
over potentially small scales (such as the formation region of CR chondrites) are 
not incorporated into the current approach. Similarly, the dust/gas ratio is also 
well known to alter the oxidat ion state of a system and the stable condensates 



(see Ebel and Grossman . 2000l and references therin). Such variations in the ox- 



idation state may prove to be necessary to produce the parent bodies of certain 
meteorite classes, such as the enstatite meteorites, which are not produced under 
the current equilibrium assumption. Enstat ites contain a very reduced miner- 
alogy, suggesting a C/0 value above solar ( Larimer and Bartholomav . 1979I ). 



Such regions and localized variations are neglected from the current approach 
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as they are not yet well understood or simulated. The migration of material 
within the disk can also act to alter the oxidation state of the solid material. For 
example, the redistribution of water within the inner 5AU of the disk over time 
via diffusion and adve ction has been found to signi f icantl y alter the oxidation 
state of the disk itself ( Pasek et al. . l2005HCvr et"aI1 . ll999l) . producing both re- 
ducing and oxidizing regions of various widths and locations (depending on the 
initial conditions). As a result, the redox state of the solid material should be 
drastically altered as the disk evolves. Until we are better able to model these 
variations, we are unable to examine their effects on the final simulated plan- 
etary compositions produced. It is also worth mentioning that the extremely 
reduced nature of the solid material may also be increased to some extent by our 
present inability to simulate the olivine and pyroxene solid solutions. This acts 
to limit the amount of Fe that can be oxidized and incorporated into silicate 
species, thus producing more reduced compositions. 

Comparison of simulated oxidation values to those of the terrestrial planets is 
difficult, either due to lack of direct information (Venus) or a varying redox state 
(Earth and Mars). Earth has a reduced m etallic core and more oxidized mantle 
(near QFM) (| Joned . |200 j : iRollinsonl . l2007t ) . Although the SNC meteorites span 
a large range of oxidation states (IW to I W-|-2.jl), the martian interior itself may 



have an oxidation state close to IW— 1 (jJonesl . 120041 ). Defining the "average" 
oxidation state of the Earth is complicated by va riations due to te ctonic settings, 
especially for the crust and uppermost mantle (|Rollinso Furthermore, 
it is believed that the redox state has varied through time (e.g. lGalimovll2005f ). 
In order fo r core f ormat ion to occur, a mantle composition of lW-2 to IW-3 
is required ( Joned . [2004h . However, the mantle must have achieved its present 
oxidation level (~QFM) relatively early in the Earths history as little temporal 
variation is o bserved i n the oxidation state of basalts younger than ~3.5 - 3.9 
billion years ()Delanol . [200lh . As a result of these issues there is currently no 
way of accurately determining an oxidation state for the primitive terrestrial 
planets (and thus the time of the 'best fit' disk conditions implied by it) at the 
stage of formation that we are simulating here. However, the generally reduced 
nature of the planets produced in the current simulations is in agreement with 
the hypothesis of an initially reduced Earth and Mars before planetary processes 
altered the redox state of their upper layers. 



3.4- Orbital properties of Jupiter and Saturn 

Although differences in the orbital and dynamical properties of the final 
planets are produced by the two different types of simulations considered here 



( O'Brien et al. . 20061 ). the differences between the CJS and EJS simulations in 



terms of both their bulk elemental abundances and their geochemical ratios for 
rock forming elements are negligible. Neglecting the small volatile rich planet 
formed in the asteroid belt in simulation EJS3, comparable bulk compositional 



■^IW refers to the iron-wustite buffer. IW+2.5 refers to an oxygen fugacity 2.5 log units 
above that of the buffer. 
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trends can be seen in both the CJS and EJS simulations. This clearly indicates 
that in the case of terrestrial planet formation within the Solar System, the bulk 
composition of the final terrestrial planets formed is not highly dependent on 
the orbital properties of Jupiter and Saturn. Difference between the two types 
of simulation emerge when we consider the delivery of hydrous species (Section 

Once again, this is a result of the fact that the majority of planet forming 
elements have fully condensed out of the nebula by 1200K, co rresponding to 



midpl ane radii within 0.75AU from the host star for the models of lHersant et al 



(|200lh and do not change greatly in their relative weight percentage values over 
the remainder of the simulation region. As a result the majority of the solid 
mass within the present simulations has solar elemental abundances. This is 
consistent with the idea that the equilibrium chemical composition of the Solar 
nebula is not highly zoned in that it does not contain radially narrow zones 
of vastly different compositions. Instead, the vast majority of the solar disk is 
dominated by material composed of pyroxene (MgSiOs), olivine (Mg2Si04) and 
metals (Fe and Ni). 

Figure [5] shows the weight percentage (in solid material) of the key planet 
building elements O, Fe, Mg and Si, in addition to the mass distribution of the 
dynamical simulations. It can be seen that the majority of mass is expected to 
have similar relative elemental abundances. The small variations observed in the 
weight percent values are caused by increases in the solid portion of disk mass 
generated by the condensation of FeS at lAU and serpentine (Mg3Si205 (011)4) 
at 3.5AU. Thus small variations in the feeding zones produced by different 
orbital properties of Jupiter and Saturn within the CJS and EJS simulations 
are unlikely to produce major differences in final bulk composition. 

3.5. Variations with Time 

Figure [6] shows the variation with time for weight percentage values of several 
key planet building elements in the planets produced by the CJSl and EJSl sim- 
ulations. It can be seen from Figure [S] that very little variation in composition 
(<5 wt%) occurs during the formation process with the bodies obtaining their 
final compositions relatively early during accretion and displaying only minor 
deviations over time. The same trend is also observed for the six other simula- 
tions not shown. This implies that the planets within these simulations formed 
homogenously (i.e. from material with similar composition to the final planet). 
If formation did occur in this manner, then the presence of a 800-lOOOkm deep 
magma ocean would be required during core formation in order to produce the 
sidero phile element abundances observed in the crust and upper mantle of the 
Earth (jPrakel . [2OO0I ). 



However, it is likely that the homogenous accretion observed here is due to 
the 'snapshot' approach we have taken when determining the composition of 
solid material within the disk. In these simulations, we have considered the 
composition determined by disk conditions at just seven discrete times. In real- 
ity, it is quite likely that the composition of solid material will change over time 
as it undergoes migration and experiences other disk and stellar processes (such 
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as the redistribution of water ( Cvr et al. . 19991 )). Such changes are not captured 



in our current approach. As such, more detailed chemical models incorporating 
temporal variations in the composition of solid material are required in order to 
test the homogenous accretion produced here. 

3.6. Late Veneer 

An alternative hypothesis to homogenous accretion is heterogeneous accre- 
tion where the composition of material accreted changes significantly with time. 
Under this hypothesis, the highly siderophile elem ent distribut ion within the 
Earth is explained via accretion of a "late veneer" (|Choul . ll978h . Referring to 



the last ~1% of mass accreted by the Earth, this material would need to contain 
chon dritic abundances o f the s iderophile elements, and essentially no metallic 



iron (jPrake and Righteij . 120021 ) . Although the CI and CM carbonaceous chon- 
drites are sufficiently oxidized and contain limited metallic Fe, their Os isotopic 
abundances are not in agreement with values required for the late veneer. On 
the other hand, H-chondrites do posses Os isotopic ratios of the correct value 
but are not sufficiently oxidized. As such, we currently have no samples in the 
meteorite record that rep resent a possible source for the late veneer material 



meteorite record tnat rep resent a possible source lor tne late veneer material 
( Drake and Righteill2002f ). However, this is not surprising, given the limited 



number of chondritic parent bodies in our meteorite collection (^^15) compared 
to the thousands of planetesimals that accreted to form the Earth. 

As expected from the previously observed homogenous accretion, the late 
veneer of the present simulations is similar in composition to the final plane- 
tary abundance. Here the late veneer is taken to be the material accreted after 
the last impact by a projectile with mass > y^embryo- This material is not 
highly oxidized and is primarily composed of olivine, pyroxene, metallic iron, 
troihte (FeS), diopside (CaMgSi206), nickel and albite (NaAlSiaOs). This com- 
po sition is similar to th at of the ordinary chondrites. Furthermore, as noted 
in O'Brien et al. I (|2006l) . the planets produced in the EJS simulations accrete 



an average of ~10% of the final planetary mass as a late veneer, an order of 
magnitude above the predicted amount. Thus it can be seen that the current 
simulations do not successfully reproduce the late veneer material. However, as 
discussed in sections 13.11 and 13.51 migration of material and temporal variations 
in the composition of solid material accreted have the potential to drastically 
alter the composition of the late veneer. Such variations are not currently cap- 
tured by the current simulations. 

3.7. Hydrous Species 

Differences between the CJS and EJS simulations emerge when considering 
the delivery of hydrous species to the final planets. In the present simulations, 
"hydrous species" refers to water ice and the aqueous alteration product ser- 
pentine (specifically clinochrysotile Mg3Si205(OH)4). Althou gh silicate melts 



on th e bodies themselves may contain water (up to 15 wt% (jHamilton et al 
Il964[ ))and thus also be an important source of hydrous material, they are cur- 
rently neglected from this study as we are only focussing on the equilibrium 
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condensation composition of the system. As stated in O'Brien et al. ( 20061 ). 



the planets formed in the EJS simulations do not accrete significant amounts of 
volatile rich material from beyond 2.5AU (with the exception of the same planet 
formed in the E JS3 simulations as previously mentioned) . Thus for disk condi- 
tions before 1x10^ years, these planets are found to contain no hydrous species. 
For disk conditions after this time, however, all terrestrial planets produced in 
the EJS simulations are found to have significant water components, although 
generally less than those of the CJS terrestrial planets. The appearance of water 
and serpentine at later times is largely due to the migration of the ice line to 
within the region from which source material is obtained. The planets formed 
via the CJS simulations, however, all contain significant amounts of serpentine 
and water ice for disk conditions at t = 5x 10^ years. At this time, seven planets 
include hydrous material (<O.OO8M0), incorporated via the accretion of one 

to four planetesimals each containing a minor amount of serpentine. Only one 
planet contained a large amount of serpentine (>O.O1M0). All eight planets 
accreted the hydrous material before the late veneer and during the stage where 
large, violent impacts were occurring. As such, it is expected that a significant 
component of this hydrous material would be vaporized during later impact 
events and subsequently lost from the final planetary body or incorporated into 
the planetary core. Only four planets (including two Earth analogs) accreted 
serpentine (0.001 - O.OOSM^) as part of their late veneer. Given the relatively 

late delivery of this material, it is believed that the majority would be retained 
during later planetary processing. As one would intuitively expect, the amount 
of hydrous material accreted by a planet increases with increasing planetary 
semimajor axis. 

Assuming that all of the hydrogen accreted as hydrous species is converted to 
water, the current CJS simulations are producing planets containing 0.6 to 24.8 
Earth ocean masses of water for disk conditions at 5x10^ year^. If we assume 
that Venus initially possessed a similar amo unt of water to the E arth while Mars 
contained 0.06 - 0.27 Earth oceans of water ( Lunine et al. . 20031 ). then it can be 



seen that for all of the terrestrial planets considered, the present simulations are 
producing planets with sufficient water to avoid the need to invoke other large- 
scale exogenous water delivery sources such as cometary impacts. Of course, 
these values should be considered to be upper limits on the amount of total 
water present within the planet as potentially large amoun ts (>50%) may be lost 
during the accretion process (jCanup and Pierazzo . |2006[ ). Additionally, some 



will undoubtedly be lost by photodissociation, Jean's escape and (definitely in 
the case of Earth) by the formation of organic species. Primordial accretion of 
water (as predicted here) would result in H being accreted to the core of the 
planet via increased H partion i ng into molten FeHx from a hydrated silicate melt 
( Okuchi and Takahashil . l l99ilOkuchil . ll99it . This would again act to decrease 



the amount of water available on the planetary surface. Although values for 
-'l Earth ocean mass of water = 1.4xl0^^kg = 2.34xlO"*M0 
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the amount of water expected to be lost by these processes are unknown, it is 
still likely that a significant portion of the initial water will be retained within 
or on the crust. The presence of water would also have a variety of effects on 
the nature of the proto-earth. It would alter the phase relations for the key 
components of a magma ocean, modify the physical properties of the melt (such 
as its comp ressibility) and also the crystal settling properties of such a melt 
( Abe et alllioooh . 



Temporal variations in the composition of the disk are likely to increase the 
amount of hydrous material accreted by the planets due to water-rich material 
being stable over a drastically larger fraction of the disk. These time-dependent 
variations can be observed in our current approach as simulations undertaken 
for conditions at 3x10^ years produce terrestrial planets containing up to 1200 
Earth ocean masses of water, well above the levels observed for simulations under 
disk conditions at 5x10^ years. The current approach also does not account 
for the possibility of water delivery yia adsorption of wa,ter onto solid grains 
later incorporated into planetesimals ( Drake and Campinsl . 2006t Stimpfl. et al 



20041) . As all of the solid material considered in the current simulations would be 
bathed in H2O vapor, it is possible that a significant amount of water could be 
delivered to the final planets via adsorption that is not presently accounted for. 
Thus it appears likely that the terrestrial planets form "wet" with a sizeable 
portion of their primordial water delivered as a natural result the planetary 
formation process. 

This conclusion does not consider the resulting D/H ratio of the accreted wa- 
ter, a key constraint on any water delivery hypothesis, as isotopic abundances 
can not be determined in our current approach. Primordial D/H ratios are 



know n to vary greatly, from 2— 3x 10 ^ for protostella r hydrogen dLecluse and Robert . 



1994 ) to 9 X 10 ^ for aqueous inclusions in meteorites ( Deloule and Robert 19951 ) 



Similar variations can be seen in solid bodies with carbonaceous chondrites con- 
taining values close to that of the Earth (1.5xl0 ~'^) while Mars appears t o 



be enriched in D with D/H values around 3x10""* (|Drake and Righteil . l2002f ). 
Thus it is difficult to make detailed predictions about the possible D/H ratios of 
the simulated planets. However zeroth order predictions can be made if we as- 
sume that the D/H ratio decreases linearly with increasing semimajor a xis from 
the o bserved Martian value to that of Jupiter's atmosphere ("2.1 xlO~^) (|Robertl . 



I2OO6I ). As all of the hydrous species in the current simulations are produced 
beyond 3.6AU (for disk conditions at t = 5x10^ years), we can assume that the 
D/H ratio will be less than 1.4x10""*. Thus it is expected that the D/H ratios 
for the simulated planets will be less than the currently observed planetary val- 
ues. However, given the large degree of processing this material is expected to 
experience both during and after planet formation, it is likely that the planetary 
D/H values will increase as H is preferentially lost from the system. 

3.8. Volatile Loss 

The loss of volatile material (pres umably through vapo r loss) in large impact 
events can potentially be significant. lO'Keefe and Ahrend (1977) found that for 



impacts between two anorthosite bodies the volume of vapor produced was less 
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than that of the projectile for impact speeds of 15 km/s and increased up to 
10 times that of the projectile at speeds of 45 km/s. This increased up to 40 
times the volume of the proj ectile for a collision betwee n an iron body and 
anorthosite body at 45 km/s ( O'Keefe and Ahrend . Il977 ). Production of this 
vapor plume can result in significant devolatization of the target body. For 
example, 2 to 60% of water is lost fr om a water-bearing target body during 
an impact (jMartv and Yokochil . l2006h while impacts into serpentine samples 
produced more than 40 % devolatization for peak shock pressures above ~15 GPa 
( Tvburczv et al. . 200l[ ).Similarlv. shock heating studies of Mur chison meteorite 



samp les have shown losses of up to 76% of H from the target (jTvburczv et al 
19951) . 



Although these losses of volatile material arc important, detailed studies of 
such a loss throughout the duration of the terrestrial planet formation process 
have not been undertaken. Affects such impacts into partially molten bodies, 
loss of elements an species other than H and water, effective retention of the 
vapor plume by the growing body and subsequent re-condensation may become 
important during planetary accretion. As such, we are limited to making only 
first-order approximations of the amount of each element lost from the final 
planet due to impacts. To do this, we need to determine both the amount of 
the final body that is molten and/or vaporized after each impact and the amount 
of each element that woul d be lost from the molten phase. Using equation 9 from 
Tonks and Melosh (|l993l) . the volume of melt produced by the initial Shockwave 
of an impact, Vm, for a hemispheric melt model is given by: 



V„,, = V, 



PP 



pro'j 



PtCOS 



1/2,; 



3/2 



(3) 



where Yproj is the volume of the projectile, pp is the density of the projectile, pt 
is the density of the target, i is the impact angle (as measured from the vertical), 
Vi is the impact speed and is the minimum impact speed needed to produce 
melting. 

Thus it can be seen that to first order: 

oc (4) 



Since the fraction of the planet that is molten (fm) is also directly propor- 
tional to the volume of the planet that is molten, we can deter mine the melt frac- 
tion p roduced by each impact by scaling the fraction given in lTonks and Melosh 
for a Vj = 15 kms ^ i mpac t with the v^, Mproj and yitarget values deter- 
mined by the lO'Brien et all ( 2006^ models (assuming a uniform density for both 
the target and the projectile). Although a 10-30% difference in the melt volume 
pr oduced by the truncated sphere and hemisphere melt models was identified 
bv lTonks and MeloshI (|l993l) . the scaling applied here is acceptable for the goals 
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of the current study. Note that this approach assumes constant values for i 
and v™-. Although these values will undoubtedly change, we are unable to pro- 
vide any limitations for t hem at this time. As such, we adopt the values from 
Tonks and MeloshI (|l993l ). 



The amount of each element assumed to be lost from the molten phase is 
based on the condensation temperature of the element. The most volatile el- 
ement present in solid form (C, primarily present as CH4.7H2O in the Solar 
nebula) was set to lose 100% of its mass from the melt. For all other elements, 
the percentage lost was determined from the observed elemental depletions in 
ordinary chondrites (as compared to C I chondrites') with respect to the 50% 
condensation temperature as taken from David ( 20061 ). CI chondrites are widely 
believed to be the most chemically primitive meteorites. Thus the elemental 
depletions relative to CI chondrites observed in other meteorite classes are be- 
lieved to be due to processing of material within the solar protoplanetary disk 
and as such serves as an excellent proxy for the loss of volatile material. The 
percentage of each element assumed to be lost in each impact event is shown 
in Table O Note that the condensation temperature of O is taken to be the 
condensation temperature of silicate, the dominant form of O throughout the 
majority of the disk. The amount of each element lost was then determined after 
each impact event and all lost volatile material was assumed to be permanently 
removed from the planet into the surrounding disk. 

This approach is obviously based on several broad assumptions and can only 
provide order of magnitude approximations for the loss of volatile material. In 
addition to assuming details about each impact event (such as cos-^/^z), this 
approach assumes that the target body has completely cooled and resolidified 
between impacts (i.e. each impact oc curs with two cold solid b odies). A hotter 
body produces a larger melt fraction (jTonks and MeloshI . Il993h and is expected 
to result in greater loss of material from each impact event. The current simu- 
lations are b ased on the equation of st ate for dunite obtained from the AENOS 
package fsee lTonks and MeloshI . Il993l for Hugoniot details). Variations in this 
may also alter the degree of melting experienced by a body. Furthermore, we 
currently only consider melting produced by the initial Shockwave as it moves 
through the body. Detailed hydrocode simulations need to be undertaken in 
order to examine how volatile losses would vary under more realistic conditions 
including differentiation (or lack thereof) of the target body, atmospheric losses, 
impacts into a still-molten embryo and recapturing of ejected material. Finally, 
the current dynamical models assume perfect accretion for each of the impact 
events. To obtain more realistic simulation conditions, loss of volatile mate- 
rial needs to be incorporated into the simulated impact events. However, this 
approach is currently computationally very demanding and is not feasible. 

As expected, loss of material through impacts reduced the amount of volatile 
species present in the final planetary bodies but did not significantly alter the 
abundance of more refractory species. This loss of volatile elements can best be 
seen in the planet normalized abundances shown in Figure [7| Clear reductions 
in the amount of volatile elements (specifically Na and S) can be seen, while neg- 
ligible changes are produced in the abundances of the more refractory elements 
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(such as Fe, Cr and Mg). These reductions in the abundance of the most volatile 
elements are unable to produce final planetary abundances in exact agreement 
with the observed planetary abundances within the Solar System but they do 
represent a substantial improvement, and indicate that impact-induced melting 
and vaporization (and the associated loss of material) is an important factor in 
determining the bulk elemental abundances of volatile species within a planet. 

The average fraction of the planet melted in each individual impact event is 
< 5% for both simulations (2.8% for the CJS simulations, 3.5% for the EJS sim- 
ulations). Only one simulation (EJSl) produced an impact event large enough 
to melt and/or vaporize the entire target body. This event was the 10th impact 
occurring on this body and thus occurred early in the dynamical formation sim- 
ulation (t = 9.72x10^ years for the dynamical simulation). As such, it is likely 
that although the entire body would have been disrupted by the event, a solid 
body may have reformed from the remaining material afterwards. More detailed 
simulations are required to determine the full effects of large-scale impacts such 
as this one on the terrestrial planet formation process. 



3.9. Solar Pollution 

Insignificant amounts of pollution of the Solar photosphere occurred during 
the terrestrial planet formation simulations. A maximum of 0.135M q^ of solid 
material was added to the Sun in the CJS simulations, while the EJS simulations 
contributed a maximum of l.llM^. The ensemble-averaged resulting solar 
abundances are shown in Table [HI The addition of solid material generated in 
the CJS simulations produced no observable enrichment in the Solar spectrum. 
The EJS simulations, however, did produce a minor enrichment of up to 0.02 
dex for Ti and Al, 0.01 dex for C, N, Na, Mg, Al, Si, P, S, Ca, Cr, Fe and 
Ni and no enrichment in O. Nonetheless, this enrichment is clearly not large 
enough to be definitively detected with current spectroscopic studies as many 
such studies produce stellar abundances with errors equivaleii t to or larger than 
the expected enrichment (e.g. ±0.03 for Fischer and Valenti ( 20051 )1. As such, 
any Solar pollution produced by terrestrial planet formation with the Solar 
System is believed to be negligible. 



4. Discussion 



The simulations successfully produce terrestrial planets that are in excellent 
agreement with the terrestrial planets of the Solar System, in terms of both their 
dynamics and their bulk elemental abundances. Although current simulations 
are unable to capture the finer details of planetary formation (such as the com- 
position, and sometimes the amount, of the late veneer material), the success 
of these simulations on a b roader scale provides us with increased confidence 



in the dynamical models of lO'Brien et al.l (l2006f). as well as other comparabl 



high-resolution N-body simulations such as iRavmond et~all (l2006l [20071 [200 
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(earlier, lower- resolution simulations such as Chambers and Wetherill ( 1998f ). 
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Chamber^ ( 2001 ). and Raymond et alj ( 2004 ) would likely also provide a rea- 



sonable fit in terms of bulk chemistry, but are unable to fully match properties 
of the terrestrial planets such as their masses and number, formation timescale, 
low dynamical excitation, etc). Furthermore, it also serves to validate the ap- 
proach utilized here to combine detailed dynamical and chemical modeling to- 
gether. This will allow for reliable application of this approa ch not only t o other 
dynamical models but also to other planetary systems fe.g. iBondl . 120081 ). 

The final bulk elemental compositions of the simulated terrestrial planets 
for rock forming elements are not strongly influenced by the orbital properties 
or evolution of Jupiter and Saturn as can be seen by the strong similarities be- 
tween the CJS and EJS simulations. The only significant differences occur in the 
amount of water rich material accreted onto the final planets and the amount 
of late veneer material delivered after the last major impact. This suggests that 
the bulk chemical evolution of the terrestrial planets is to a large extent inde- 
pendent of the evolution of the giant planets. This conclusion, however, is only 
valid for late stage and in situ formation after the giant planets have formed 
and undertaken the vast majority of their migration. Simulations are currently 
running to examine the effects on planetary composition of formation occur- 
ring during giant planet migration. As Jupiter and Saturn are not believed to 



have undergone extensive migration (jGomes et al.l . 120051 : iLevison et al.l . 12005 



Morbidelh et al.U2005l ). the issue of migration is not a significant one for the 



Solar System. However, it will be an important issue for extrasolar planetary 
systems where many are thought to have experienced large amounts of migra- 
tion. 

The delivery of water and other volatile species to the final terrestrial planets 
appears to be a normal outcome of terrestrial planet formation within the CJS 
simulations, implying that it is normal for such planets to accrete significant 
amounts of water. As such, it negates the need for large-scale delivery of water 
and other volatile species by such exotic processes as cometary impacts. The 
EJS simulations, on the other hand, only produce terrestrial planets with a 
significant amount of water for disk conditions at times of 1x10^ years and 
greater. Thus in order to produce terrestrial planets with water accreted during 
the formation process within the EJS configuration, time-varying compositions 
need to be included within a single simulation, allowing for an increased hyrdous 
material content at later times. It should also be noted that the terrestrial 
planets produced in the CJS simulations are likely to always accrete more water 
than their EJS companions. Terrestrial planets within the CJS simulations 
accrete more material from the outer, more water enriched regions of the asteroid 
belt, thus making them inherently water-rich compared to the terrestrial planets 
produced in the EJS simulations. 

Later planetary processing (especially within the mantle) would have re- 
sulted in dissociation of some of the primordial water, causing a significant 
amount of H to migrate to the core of the Earth or alternatively to be lost from 
the planet. This migration, however, would have left a large amount of OH in 
the mantle, thus increasing its redox state. This process allows us to explain 
the evolution of the Earth's redox state over time from the initially reduced 



23 



state produced by the current simulations to the present stratified redox state. 
This conclusion, however, requires a more detailed experimental understanding 
of the efficiency of H partioning within the mantle and core of the Earth, along 
with the evolution of mantle processes and mixing. 

It is interesting to note the order of magnitude difference between the CJS 
and EJS simulations in the amount of solid material added to the Sun dur- 
ing terrestrial planet formation. This difference in accreted material is due to 
stronger resonances produced when Jupiter and Saturn are in more eccentric 
orbits. These resonances can quickly increase the eccentricity of a body to 
effectively driving them into the Sun and resulting in an increased amount of ac- 
creted material. This may potentially have great bearing on extrasolar planetary 
host stars as many known extrasolar planets are currently in eccentric orbits. If 
such an orbit always results in a greater amount of solid material being accreted 
by the host star, then it may give more weight to the pollution hypothesis which 
has previously been suggested to explain the observed [ Fe/H] enrichment (e.g. 



lLaughlin|[2000l : I Gonzalez et al.l 200lUMurrav et al.ll200lh . However, no correla 



tion has been found between known planetary eccentricity and the metallicity 
of the host star (e g. iReidl . I2OO2I : ISantos et all . i2003aHFischer and Valentil . 12005 ; 
Bond et al. ,[2Q0i), suggesting that this is not a strong effect. We determined 



that a negligible amount of solid material was added to the Solar photosphere 
during terrestrial planet formation and that no observable elemental enrich- 
ment would be produced. While similar simulations are needed for extraso- 
lar planetary systems, the current simulations support the conclusion that the 
observed metal enrichment in extrasolar host stars is primordial in origin, es- 
tablished in the giant molecular cloud fro m which these sys t ems fo r med as has 



been concluded by other studies (such as ISantos et al.ll200lL l2003al lbl. l2005l and 



Fischer and Valenti 2005f) . Of course, this does not rule out the possibility that 



extrasolar planetary host stars may have in fact accreted a giant planet during 
the formation and migration process and thus display higher levels of pollution. 

Finally, biologically important elements are obviously of great interest, espe- 
cially for the Earth. Of the six major biogenic elements (H, C, N, O, S and P), 
four are accreted in excess by the planets during their formation (H, O, S and 
P). Only C and N are not accreted during planetary formation. Both species 
are primarily present (under the assumption of equilibrium) as solids within the 
Solar nebula as clathrates (methane or ammonia trapped in a water ice lat- 
tice) and organics. As these species only form in the outermost regions of the 
disk considered in the current simulation (beyond ~3.6AU for disk conditions 
at 5x10^ years), they could only be delivered to the final planets via cometary 
impacts, migration or temporal variations within the disk or some combination 
of all three. Therefore delivery of material from the outer regions of the disk 
are necessary in the current models for life to be able to develop. However, it 
is entirely possible that both C and N could be present at low levels in solid 
solutions and/or liquids within the actual Solar System. Furthermore, both the 
CM and CR chondrites arc known to contain up to 5 wt.% organic C. Such 
compositions are not captured in our current assumption of equilibrium driven 
condensation, indicating the presence of another mechanism, acting to incor- 
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porate this material into carbonaceous asteroids. As such, small amounts of C 
and N may be delivered during the normal accretionary process, thus reducing 
the requirement for external delivery of substantial amounts of material. 

We did not perform a comparable chemical analysis for other recent types 
of terrestrial planet formation simulations, such as N-body simulations incorpo- 
rating the e ffects of tidal drag and sweeping secular resonances during gas disk 
dissipation (iThommes et al . 2008), or the hybrid N-body /co agulation code of 
Bromlev and Kenvon ( 20061 ) and Kenvon and Bromley ( 20061 ). However, we ex- 



pect that such simulations would also provide a reasonable fit to the terrestrial 
planets in terms of bulk chemistry. As discussed in Sec. [31 the equilibrium 
chemical composition of the Solar nebula is not highly zoned, and the majority 
of the solid mass within the disk (at least outside of ~0.75 AU) has elemental 
abundances that are basically solar. While pro cesses that are treated in these 
other models [eg. secula r resonance sweeping in^Tho mmes et all (2008h and gas 



drag o n planetesimals in lBromlev and Kenvo n (2006) and Kenvon and Bromlev 



( 2006f) ] may lead to a somewhat increased degree of radial migration of material 



during terrestrial planet formation, that material will still have a fundamentally 
similar composition throughout the disk, such that the final bulk planet chem- 
istry is unlikely to be significantly different than what we find in this work. More 
subtle changes may occur, for example with regards to the amount of hydrated 
material delivered, or the amount and chemistry of the late veneer material, 
which would be interesting to explore in future work. 

5. Summary 

Bulk eleme ntal abundances have been determined for the simulated terres- 
trial planets of lO'Brien et al. I (|2006h . These abundances are in e xcellent agree- 



ment w ith observed planetary values, indicating that the models of lO'Brien et al.l 



( 2006f ) are successfully producing planets comparable to those of the Solar Sys 



tem in terms of both their dynamical and chemical properties, adding greater 
weight to their predictive properties. Simulated redox states are also in agree- 
ment with those predicted for the early Earth. Although differences do exist 
between the observed and predicted geochemical ratios, these are believed to be 
a result of our assumption of equilibrium controlled compositions. Additionally, 
the current simulations are unable to successfully reproduce the accretion of a 
late veneer of material by the early Earth, in terms of the chemistry and, in the 
case of the EJS simulations, also the amount. 

Significant amounts of water are accreted in the present CJS simulations, 
implying that the terrestrial planets in these simulations form "wet" and do 
not need significant water delivery from other sources. Terrestrial planets pro- 
duced in the EJS simulations only accrete significant amounts of water for disk 
conditions at 1x10^ years and beyond. N and C, however, do still need to be 
delivered to an early Earth by some other process in order for life to develop. 

Additionally, the bulk abundances of the major rock-forming elements for 
the final planets in the current simulations are not strongly dependent on the 
orbital properties of the giant planets, with the CJS and EJS simulations both 
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producing comparable results. This suggests that although the orbits of Jupiter 
and Saturn are clearly of great dynamical importance to the evolution of the 
terrestrial planets, and can have a significant influence in terms of water delivery, 
they do not exert such a large influence over the chemical evolution of the 
terrestrial planets with regards to the rock-forming elements. 

Finally, the pollution of the outer layers of the Sun via solid material during 
planetary formation produces a negligible photosphcric elemental enrichment. 
Assuming similar levels of pollution in other planetary systems, this in turn 
implies that the high metallicity trend observed in extrasolar planetary systems 
is in fact primordial. 
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T able 1: Propert i es of s imulated terrestrial planets produced in the Solar System simulations 
of lO'Brien et al.l l|2006l V Numbering starts at 4 and increases with increasing distance from 
the Sun. CJS denotes the results of the circular Jupiter and Saturn simulations while EJS 
indicates the results of the eccentric Jupiter and Saturn simulations. 
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Table 2: HSC Chemistry input values for Solar System Simulations. All values are in moles. 
All species are assumed to initially be in their elemental and gaseous form with no other 
species present within the system. 



Element 


Abundnace 


H 


1.00 xlO^^ 


He 


8.51 xlQio 


C 


2.45 xIQS 


N 


6.03 xlO^ 





4.57 xIQS 


Na 


1.48 xIQS 


Mg 


3.39 xlO^ 


Al 


2.34 xIQS 


Si 


3.24 xlO^ 


P 


2.29 xlO^ 


S 


1.38 xlO^ 


Ca 


2.04 xlO^ 


Ti 


7.94 xlO^ 


Cr 


4.37 xlO^ 


Fe 


2.82 xlO^ 


Ni 


1.70 xlO*^ 
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Table 3: T^ox condensation foi' ^ g^^ with Solar photosphere composition and at P = 10~^ 
bar. Ini t ial ph ase for each element is also listed. Values are provided from the models of 
iLodder j ll2003h for comparison. 



Element T5o% condensation (K) Initial Phase 





This Study 


Lodders (2003) 


Species 


Formula 


Al 


1639 


1665 


Hibonite 


CaAliaOig 
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CH4.7H2O 
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FeS 


Si 


1329 
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Gehlenite 


CaaAlaSiOr 


Ti 


1580 


1573 


Perovskite 


CaTiOa 
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Table 4: Predicted bulk planetary abundances for the terrestrial planets of the lO'Brien et al simulations for 

disk c onditions at t = SxlO^yr. All values are wt% of the final planet. CJS denotes the simulations of lO'Brien et all 
l|2006f ) with Jupiter and Saturn in the circular orbits predicted by the Nice model while EJS denotes the results of 
the simulations with Jupiter and Saturn in their current elliptical orbits. Planetary numbers start at 4 and increase 
with increasing distance from the Sun. 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




t=5xl0^ years 


CJSl-4 


0.01 


14.68 


32.20 


2.90 


27.73 


1.37 


1.99 


0.42 


1.76 


0.39 


0.12 


0.09 


16.34 


0.00 


0.00 


0.90 


CJSl-5 


0.01 


14.86 


31.36 


2.97 


28.42 


1.19 


1.74 


0.60 


1.80 


0.41 


0.14 


0.08 


16.44 


0.00 


0.00 


0.90 


CJSl-6 


0.01 


14.09 


30.86 


7.03 


26.94 


1.08 


1.57 


0.58 


1.71 


0.39 


0.13 


0.08 


15.56 


0.00 


0.00 


0.91 


CJS2-4 


0.00 


15.06 


32.28 


1.11 


28.07 


1.65 


2.39 


0.20 


1.79 


0.39 


0.09 


0.11 


16.88 


0.00 


0.00 


0.89 


CJS2-5 


0.00 


14.74 


31.06 


3.97 


28.19 


1.13 


1.65 


0.59 


1.79 


0.41 


0.14 


0.08 


16.28 


0.00 


0.00 


0.91 


CJS2-6 


0.00 


14.61 


31.53 


4.07 


27.90 


1.14 


1.66 


0.57 


1.77 


0.40 


0.13 


0.08 


16.15 


0.00 


0.00 


0.90 


CJS2-7 


0.08 


13.83 


31.71 


7.28 


26.42 


1.06 


1.54 


0.57 


1.67 


0.38 


0.13 


0.07 


15.27 


0.00 


0.00 


0.91 


CJS3-4 


0.00 


14.87 


32.26 


2.12 


27.92 


1.45 


2.11 


0.34 


1.78 


0.39 


0.11 


0.10 


16.57 


0.00 


0.00 


0.90 


CJS3-5 


0.00 


14.48 


31.21 


4.96 


27.70 


1.11 


1.62 


0.59 


1.75 


0.40 


0.13 


0.07 


15.99 


0.00 


0.00 


0.91 


CJS3-6 


0.12 


13.70 


32.35 


7.20 


26.16 


1.05 


1.53 


0.56 


1.66 


0.38 


0.13 


0.07 


15.12 


0.00 


0.00 


0.91 


CJS4-4 


0.00 


14.86 


31.72 


2.61 


28.18 


1.34 


1.95 


0.44 


1.79 


0.40 


0.12 


0.09 


16.51 


0.00 


0.00 


0.90 


CJS4-5 


0.05 


14.01 


31.47 


6.73 


26.78 


1.07 


1.56 


0.58 


1.70 


0.39 


0.13 


0.07 


15.47 


0.00 


0.00 


0.91 


EJSl-4 


0.00 


15.11 


31.73 


1.34 


28.82 


1.31 


1.90 


0.59 


1.82 


0.42 


0.14 


0.09 


16.76 


0.00 


0.00 


0.90 


EJSl-5 


0.00 


14.84 


31.34 


3.40 


28.22 


1.21 


1.77 


0.41 


1.79 


0.40 


0.12 


0.08 


16.43 


0.00 


0.00 


0.90 


EJSl-6 


0.00 


14.13 


31.27 


6.47 


27.00 


1.08 


1.58 


0.58 


1.71 


0.39 


0.13 


0.08 


15.61 


0.00 


0.00 


0.91 


EJS2-4 


0.00 


14.64 


31.61 


2.84 


27.48 


1.75 


2.55 


0.33 


1.75 


0.38 


0.09 


0.12 


16.49 


0.00 


0.00 


0.89 


EJS2-5 


0.00 


15.03 


31.62 


2.12 


28.56 


1.27 


1.85 


0.46 


1.81 


0.41 


0.13 


0.09 


16.66 


0.00 


0.00 


0.90 


EJS2-6 


0.00 


14.70 


30.66 


4.55 


28.12 


1.13 


1.64 


0.60 


1.78 


0.41 


0.14 


0.08 


16.23 


0.00 


0.00 


0.91 


EJS2-7 


0.00 


13.72 


32.22 


7.36 


26.19 


1.05 


1.53 


0.57 


1.66 


0.38 


0.13 


0.07 


15.14 


0.00 


0.00 


0.91 


EJS3-4 


0.00 


15.08 


32.01 


1.11 


28.36 


1.55 


2.24 


0.40 


1.80 


0.40 


0.11 


0.10 


16.85 


0.00 


0.00 


0.90 



Continued on next page 



Tabic 4 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




EJS3-5 


0.00 


15.03 


31.30 


2.45 


28.74 


1.15 


1.68 


0.61 


1.82 


0.41 


0.14 


0.08 


16.61 


0.00 


0.00 


0.91 


EJS3-6 


0.00 


14.43 


30.65 


5.79 


27.61 


1.11 


1.61 


0.52 


1.75 


0.40 


0.13 


0.08 


15.94 


0.00 


0.00 


0.91 


EJS3-7 


1.02 


15.25 


40.93 


0.95 


18.26 


2.44 


1.57 


0.62 


0.04 


0.19 


0.09 


3.60 


14.72 


0.00 


0.35 


1.04 


EJS4-4 


0.00 


15.20 


31.99 


0.91 


28.66 


1.45 


2.10 


0.35 


1.82 


0.41 


0.11 


0.10 


16.93 


0.00 


0.00 


0.90 


EJS4-5 


0.00 


15.09 


31.94 


1.01 


28.58 


1.45 


2.10 


0.57 


1.81 


0.41 


0.13 


0.10 


16.81 


0.00 


0.00 


0.90 


EJS4-6 


0.00 


14.36 


30.02 


6.64 


27.49 


1.10 


1.60 


0.59 


1.74 


0.40 


0.13 


0.08 


15.87 


0.00 


0.00 


0.91 


EJS4-7 


0.00 


14.31 


31.04 


5.89 


27.36 


1.10 


1.60 


0.58 


1.73 


0.40 


0.13 


0.08 


15.80 


0.00 


0.00 


0.91 



Table 5: The percentage of each element assumed to be lost from the melt produced by 
individual impact eve nts. V alues are based on depletions observed in meteorite samples and 
are taken from Davi 
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% Lost 
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100.00 


N 


131 


97.69 


H 


182 


92.90 


S 


704 
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Table 6: Mean change in solar photospheric abundances produced by pollution via accretion 
of solid material during terrestrial planet format ion. All solid m aterial migrating to within 
O.IAU from the Sun during the simulations of O'Brien et al] 1I2OO6D is assumed to be accreted. 

Change is defined as AbundanCCaftcr planet formation - Abundancebeforc planet formation- 
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Final Planetary Systems 
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Figure 1: Schematic of the results of the simulations of lO'Brien et ahl l|2006l '). The horizon- 
tal lines indicate the variation between aphelion and perihelion. The vertical lines indicate 
variation in distance from the midplane due to the inclination of the planet. Numerical val- 
ues represent the mass of the planet in Earth masses. The Solar System (SS) is shown for 
comparison. 
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Figure 2: Normalized abundances for CJSl and EJSl simulated terrestrial planets. The ter- 
restrial planet each simulation is normalized to is shown in parentheses. Values are shown for 
each of seven time steps considered with the following color scheme: black = 2.5x10^ years, 
red = 5x10^ years, green = 1x10^ years, pink = 1.5x10® years, light blue = 2x10® years, 
yellow = 2.5x10® years and dark blue = 3x10® years. Left: CJSl terrestrial planets. Right: 
EJSl terrestrial planets. Referenc e Solar System planetar y abundanc es were taken from 
iMorean and Anders! ||1980| ') ( Venus') . [Kareel and Lewis! jiggsTl CEarth) and lLodders and FeelevI 
l|l997h (Marsi! 
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Figure 3: Cosmochemical ratios for all simulated terrestrial planets. Top: Al/Si v. Mg/Si for 
all simulated terrestrial planets. Bottom: Ca/Si v. Mg/Si for all simulated terrestrial planets. 
Filled circles indicate values for disk conditions at t = 2.5x10^ years while open circles indicate 
values for disk conditions at t = 5x10^ years. Values at all other times are concentrated at 
the 5x10^ years va lues and are n ot shown for clarity. Earth values are shown by filled squares 
and are taken from lKargel antTX cwis (1993) and McDonough and Sun (1995). Martian values 
are shown by a filled triangle and are taken from Loddcrs and Fcgley (1997). Venus values 
are shown by a filled diamond and are taken from iMorgan and Anders (.1980, ). 
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Figure 4: Urey-Craig oxidation state plot for CJSl and EJSl simulated planetary abundances. 
Squares indicate CJSl values. Circles indicate EJSl values. Values are shown for each of seven 
time steps considered with the following color scheme: black = 2.5x10^ years, red = 5x10^ 
years, green = 1x10® years, pink = 1.5x10® years, light blue = 2x10® years, yellow = 2.5x10® 
years and dark blue = 3x10® years. Approximate regions for several meteorite groups are 
also shown. 
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Figure 5: Distribution of solid mass and its relative composition within the Solar disk at 5 X 10^ 
years. Top: Composition (in wt%) for the solid material within the disk for O, Fe, Ms, Si and 
Al. B ottom: Initial distribution of mass within the dynamical simulations of lO'Brien et al] 
||2006|) . 
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Figure 6: Temporal variation in the elemental abundances of the final terrestrial planets 
produced by the CJSl and EJSl simulations. Variations in compositions are due to the 
accretion of embryos and planetesimals throughout the dynamical simulations. Left: CJSl 
simulation results. EJSl: EJSl simulation results. 
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Figure 7; Normalized abundances for CJSl and EJSl simulated planets showing abundances 
after volatile loss was considered. The solid line indicates the normalized abundances before 
volatile loss during impacts was considered while the dashed line indicates the normalized 
abundance once volatile loss during impacts has been incorporated. All abundances were 
determined for disk conditions at t = 5x10^ years. The terrestrial planet each simulation 
is normalized to is shown in parentheses. Reference Solar System planetary abundances 
were taken from [Morgan and Anderi lll980f) (Venus), iKargel and Lewis! ((1991) (Earth) and 
iLodders and FeelevI ||1997| ') (Mars). Left: CJSl terrestrial planets. Right: EJSl terrestrial 
planets. 
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A. Online Material 



This material is to appear online only as supplemental material. Table 1 lists 
the species (both gaseous and solid) included in the HSC chemistry calculations. 
The complete assemblage of all predicted bulk elemental abundances for each of 
the simulated planets and for each of the seven different sets of disk conditions 
examined are provided in Table [51 Values are provided as bulk wt% of the final 
planet for each set of disk conditions. Note that planetary numbers start at 4 
and increase with increasing distance from the Sun. 

Figures [8] - [12] sh ow the no rmalized abundances for each planet produced 
by the simulations of lO'Brien et al. (200(x) for simulations CJS2- 4 and EJS2 
- 4. Identical plots were included in the text for CJSl and EJSl. Normalized 
abundances are shown for each of the seven sets of disk conditions examined. 
The terrestrial planet to which each simulated planet was normalized was de- 
termined based on the semi-major axis of the terrestrial planet and is shown in 
parentheses in the upper left of each plot. 
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Table 7: Chemical species included in the equilibrium calculations of HSC Chemistry. 



Gaseous Species 



Al 


CrO 


MgOH 


PN 


AlH 


CrOH 


MeS 


PO 


AlO 


CrS 


N 


PS 


AI2O 


Fe 


N2 


s 


AlOH 


FoH 


NH3 


S2 


AIS 


FeO 


NO 


SN 




FeOH 


NS 


SO 


CH4 


FeS 


Na 


S02 


CN 


H 


Na2 


Si 


CO 


Ha 


NaH 


SiC 


C02 


HCN 


NaO 


SiH 


CP 


HCO 


NaOH 


SiN 


cs 


H2O 


Ni 


SiO 


Ca 


HPO 


NiH 


SiP 


CaH 


HS 


NiO 


SiP2 


CaO 


H2S 


NiOH 


SiS 


CaOH 


He 


NiS 


Ti 


CaS 


Mg 





TiN 


Cr 


MgH 


O2 


TiO 


CrH 


MgN 


P 


Ti02 


CrN 


MgO 


PH 


TiS 



Solid Species 



A1203 


FeSiOa 


CaAl2Si2 08 


C 


CaAli20ig 


FegP 


NaAlSisOg 


SiC 


Ti2 03 


FeaC 


Cr2Fe04 


TiC 


CaTiOs 


Fe 


Ca3(P04)2 


TiN 


Ca2Al2Si07 


Ni 


FeS 


AIN 


MgAl2 04 


P 


Fe304 


CaS 


Mg2Si04 


Si 


Mg3Si20B(OH)4 


MgS 


MgSiOa 


Cr 


H2O 




Fe2Si04 


CaMgSi2 06 
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Table 8: Predicted bulk pl anetary abundances f or the terrestrial planets of the lO'Brien et aP l|2006l ) simulations. All values are wt% of the final planet. CJS 
denotes the simulations of lO'Brien et all l l2006t) with Jupiter and Saturn in the circular orbits predicted by the Nice model while EJS denotes the results of 
the simulations with Jupiter and Saturn in their current elliptical orbits. Planetary numbers start at 4 and increase with increasing distance from the Sun. 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




t=2.5xl0'' years 


CJSl-4 


0.00 


14.84 


32.05 


2.16 


27.64 


1.65 


2.39 


0.31 


1.79 


0.38 


0.09 


0.09 


16.61 


0.00 


0.00 


0.89 


CJSl-5 


0.00 


14.86 


31.94 


1.61 


28.12 


1.66 


2.41 


0.35 


1.79 


0.40 


0.11 


0.09 


16.66 


0.00 


0.00 


0.89 


CJSl-6 


0.00 


14.73 


31.43 


3.34 


28.08 


1.27 


1.84 


0.58 


1.78 


0.41 


0.13 


0.07 


16.34 


0.00 


0.00 


0.90 


CJS2-4 


0.00 


8.29 


38.76 


0.58 


9.23 


11.24 


16.12 


0.09 


0.60 


0.12 


0.03 


0.66 


14.27 


0.00 


0.00 


0.58 


CJS2-5 


0.00 


14.89 


32.06 


0.78 


28.41 


1.77 


2.58 


0.33 


1.80 


0.41 


0.09 


0.10 


16.77 


0.00 


0.00 


0.89 


CJS2-6 


0.00 


14.60 


32.19 


2.18 


27.09 


1.94 


2.81 


0.35 


1.74 


0.37 


0.10 


0.11 


16.52 


0.00 


0.00 


0.88 


CJS2-7 


0.00 


14.57 


30.84 


4.98 


27.83 


1.12 


1.63 


0.60 


1.76 


0.40 


0.14 


0.06 


16.09 


0.00 


0.00 


0.91 


CJS3-4 


0.00 


10.61 


36.45 


1.56 


15.66 


7.70 


11.07 


0.17 


1.04 


0.20 


0.05 


0.45 


15.04 


0.00 


0.00 


0.71 


CJS3-5 


0.00 


14.88 


31.66 


2.22 


28.34 


1.43 


2.07 


0.42 


1.80 


0.41 


0.11 


0.08 


16.59 


0.00 


0.00 


0.90 


CJS3-6 


0.00 


14.24 


30.99 


5.98 


27.22 


1.18 


1.71 


0.58 


1.73 


0.39 


0.13 


0.06 


15.78 


0.00 


0.00 


0.90 


CJS4-4 


0.00 


12.41 


34.61 


1.13 


20.76 


5.35 


7.70 


0.21 


1.34 


0.28 


0.07 


0.31 


15.83 


0.00 


0.00 


0.78 


CJS4-5 


0.00 


14.69 


31.03 


4.07 


28.07 


1.19 


1.73 


0.57 


1.78 


0.40 


0.13 


0.06 


16.26 


0.00 


0.00 


0.90 


EJSl-4 


0.00 


14.48 


32.91 


0.31 


26.53 


2.69 


3.90 


0.11 


1.73 


0.35 


0.07 


0.15 


16.77 


0.00 


0.00 


0.86 


EJSl-5 


0.00 


11.97 


34.92 


1.22 


19.64 


5.90 


8.47 


0.29 


1.25 


0.28 


0.08 


0.34 


15.64 


0.00 


0.00 


0.77 


EJSl-6 


0.00 


14.79 


30.70 


4.18 


28.25 


1.13 


1.65 


0.57 


1.79 


0.41 


0.13 


0.06 


16.34 


0.00 


0.00 


0.91 


EJS2-4 


0.00 


8.74 


35.84 


0.65 


15.95 


9.42 


13.59 


0.24 


1.02 


0.22 


0.07 


0.55 


13.71 


0.00 


0.00 


0.64 


EJS2-5 


0.00 


13.57 


34.63 


0.82 


21.61 


4.32 


6.20 


0.18 


1.42 


0.28 


0.05 


0.25 


16.67 


0.00 


0.00 


0.81 


EJS2-6 


0.00 


15.20 


32.07 


0.36 


28.91 


1.49 


2.17 


0.41 


1.83 


0.42 


0.12 


0.08 


16.95 


0.00 


0.00 


0.90 


EJS2-7 


0.00 


14.19 


29.84 


7.57 


27.16 


1.09 


1.58 


0.59 


1.72 


0.39 


0.13 


0.06 


15.68 


0.00 


0.00 


0.91 


EJS3-4 


0.00 


10.60 


35.36 


0.25 


18.54 


7.51 


10.85 


0.10 


1.22 


0.23 


0.04 


0.44 


14.87 


0.00 


0.00 


0.71 


EJS3-5 


0.00 


15.27 


32.20 


0.22 


28.78 


1.56 


2.28 


0.21 


1.83 


0.41 


0.10 


0.08 


17.05 


0.00 


0.00 


0.90 


EJS3-6 


0.00 


14.60 


32.95 


1.81 


25.26 


2.41 


3.44 


0.50 


1.61 


0.36 


0.12 


0.14 


16.81 


0.00 


0.00 


0.87 


EJS3-7 


0.00 


13.79 


32.30 


6.89 


25.85 


1.26 


1.82 


0.54 


1.64 


0.37 


0.13 


0.07 


15.35 


0.00 


0.00 


0.90 
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Table 8 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




EJS4-4 


0.00 


10.86 


36.60 


0.14 


16.04 


7.84 


11.27 


0.08 


1.07 


0.20 


0.03 


0.46 


15.41 


0.00 


0.00 


0.71 


E,IS4-5 


0.00 


14.3.3 


32.71 


0.13 


26.97 


2.80 


4.05 


0.09 


1.72 


0.38 


0.04 


0.16 


16.62 


0.00 


0.00 


0.86 


EJS4-6 


0.00 


15.36 


31.87 


0.48 


29.36 


1.18 


1.72 


0.58 


1.86 


0.42 


0.14 


0.06 


16.97 


0.00 


0.00 


0.91 


EJS4-7 


0.00 


14.91 


31.31 


2.72 


28.29 


1.34 


1.95 


0.50 


1.79 


0.41 


0.13 


0.07 


16.56 


0.00 


0.00 


0.90 


t=5xl0''' years 


CJSl-4 


0.01 


14.68 


32.20 


2.90 


27.73 


1.37 


1.99 


0.42 


1.76 


0.39 


0.12 


0.09 


16.34 


0.00 


0.00 


0.90 


CJSl-5 


0.01 


14.86 


31.36 


2.97 


28.42 


1.19 


1.74 


0.60 


1.80 


0.41 


0.14 


0.08 


16.44 


0.00 


0.00 


0.90 


CJSl-6 


0.01 


14.09 


30.86 


7.03 


26.94 


1.08 


1.57 


0.58 


1.71 


0.39 


0.13 


0.08 


15.56 


0.00 


0.00 


0.91 


CJS2-4 


0.00 


15.06 


32.28 


1.11 


28.07 


1.65 


2.39 


0.20 


1.79 


0.39 


0.09 


0.11 


16.88 


0.00 


0.00 


0.89 


CJS2-5 


0.00 


14.74 


31.06 


3.97 


28.19 


1.13 


1.65 


0.59 


1.79 


0.41 


0.14 


0.08 


16.28 


0.00 


0.00 


0.91 


CJS2-6 


0.00 


14.61 


31.53 


4.07 


27.90 


1.14 


1.66 


0.57 


1.77 


0.40 


0.13 


0.08 


16.15 


0.00 


0.00 


0.90 


CJS2-7 


0.08 


13.83 


31.71 


7.28 


26.42 


1.06 


1.54 


0.57 


1.67 


0.38 


0.13 


0.07 


15.27 


0.00 


0.00 


0.91 


CJS3-4 


0.00 


14.87 


32.26 


2.12 


27.92 


1.45 


2.11 


0.34 


1.78 


0.39 


0.11 


0.10 


16.57 


0.00 


0.00 


0.90 


CJS3-5 


0.00 


14.48 


31.21 


4.96 


27.70 


1.11 


1.62 


0.59 


1.75 


0.40 


0.13 


0.07 


15.99 


0.00 


0.00 


0.91 


CJS3-6 


0.12 


13.70 


32.35 


7.20 


26.16 


1.05 


1.53 


0.56 


1.66 


0.38 


0.13 


0.07 


15.12 


0.00 


0.00 


0.91 


CJS4-4 


0.00 


14.86 


31.72 


2.61 


28.18 


1.34 


1.95 


0.44 


1.79 


0.40 


0.12 


0.09 


16.51 


0.00 


0.00 


0.90 


CJS4-5 


0.05 


14.01 


31.47 


6.73 


26.78 


1.07 


1.56 


0.58 


1.70 


0.39 


0.13 


0.07 


15.47 


0.00 


0.00 


0.91 


EJSl-4 


0.00 


15.11 


31.73 


1.34 


28.82 


1.31 


1.90 


0.59 


1.82 


0.42 


0.14 


0.09 


16.76 


0.00 


0.00 


0.90 


EJSl-5 


0.00 


14.84 


31.34 


3.40 


28.22 


1.21 


1.77 


0.41 


1.79 


0.40 


0.12 


0.08 


16.43 


0.00 


0.00 


0.90 


EJSl-6 


0.00 


14.13 


31.27 


6.47 


27.00 


1.08 


1.58 


0.58 


1.71 


0..39 


0.13 


0.08 


15.61 


0.00 


0.00 


0.91 


EJS2-4 


0.00 


14.64 


31.61 


2.84 


27.48 


1.75 


2.55 


0.33 


1.75 


0..38 


0.09 


0.12 


16.49 


0.00 


0.00 


0.89 


EJS2-5 


0.00 


15.03 


31.62 


2.12 


28.56 


1.27 


1.85 


0.46 


1.81 


0.41 


0.13 


0.09 


16.66 


0.00 


0.00 


0.90 


EJS2-6 


0.00 


14.70 


30.66 


4.55 


28.12 


1.13 


1.64 


0.60 


1.78 


0.41 


0.14 


0.08 


16.23 


0.00 


0.00 


0.91 


EJS2-7 


0.00 


13.72 


32.22 


7.36 


26.19 


1.05 


1.53 


0.57 


1.66 


0.38 


0.13 


0.07 


15.14 


0.00 


0.00 


0.91 


EJS3-4 


0.00 


15.08 


32.01 


1.11 


28.36 


1.55 


2.24 


0.40 


1.80 


0.40 


0.11 


0.10 


16.85 


0.00 


0.00 


0.90 
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Table 8 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




EJS3-5 


0.00 


15.03 


31.30 


2.45 


28.74 


1.15 


1.68 


0.61 


1.82 


0.41 


0.14 


0.08 


16.61 


0.00 


0.00 


0.91 


EJS3-6 


0.00 


14.43 


30.65 


5.79 


27.61 


1.11 


1.61 


0.52 


1.75 


0.40 


0.13 


0.08 


15.94 


0.00 


0.00 


0.91 


EJS3-7 


1.02 


15.25 


40.93 


0.95 


18.26 


2.44 


1.57 


0.62 


0.04 


0.19 


0.09 


3.60 


14.72 


0.00 


0.35 


1.04 


EJS4-4 


0.00 


15.20 


31.99 


0.91 


28.66 


1.45 


2.10 


0.35 


1.82 


0.41 


0.11 


0.10 


16.93 


0.00 


0.00 


0.90 


EJS4-5 


0.00 


15.09 


31.94 


1.01 


28.58 


1.45 


2.10 


0.57 


1.81 


0.41 


0.13 


0.10 


16.81 


0.00 


0.00 


0.90 


EJS4-6 


0.00 


14.36 


30.02 


6.64 


27.49 


1.10 


1.60 


0.59 


1.74 


0.40 


0.13 


0.08 


15.87 


0.00 


0.00 


0.91 


EJS4-7 


0.00 


14.31 


31.04 


5.89 


27.36 


1.10 


1.60 


0.58 


1.73 


0.40 


0.13 


0.08 


15.80 


0.00 


0.00 


0.91 


t=lxlO'' years 


CJSl-4 


0.63 


13.58 


34.66 


4.88 


25.95 


1.04 


1.52 


0.55 


1.64 


0.37 


0.13 


0.07 


15.00 


0.00 


0.00 


0.91 


CJSl-5 


0.19 


13.76 


32.15 


7.00 


26.31 


1.05 


1.54 


0.57 


1.67 


0.38 


0.13 


0.07 


15.20 


0.00 


0.00 


0.91 


CJSl-6 


0.83 


12.70 


36.44 


6.76 


24.26 


0.97 


1.42 


0.52 


1.54 


0.35 


0.12 


0.07 


14.03 


0.00 


0.00 


0.91 


CJS2-4 


0.07 


14.47 


31.13 


5.01 


27.70 


1.11 


1.62 


0.58 


1.75 


0.40 


0.13 


0.08 


15.98 


0.00 


0.00 


0.91 


CJS2-5 


0.12 


13.87 


31.77 


7.02 


26.48 


1.06 


1.55 


0.57 


1.68 


0.38 


0.13 


0.07 


15.32 


0.00 


0.00 


0.91 


CJS2-6 


0.56 


13.35 


34.08 


6.51 


25.52 


1.02 


1.49 


0.55 


1.62 


0.37 


0.13 


0.07 


14.75 


0.00 


0.00 


0.91 


CJS2-7 


0.76 


12.74 


36.30 


6.83 


24.33 


0.98 


1.42 


0.53 


1.54 


0.35 


0.12 


0.07 


14.06 


0.00 


0.00 


0.91 


CJS3-4 


0.75 


13.66 


35.25 


3.82 


26.10 


1.05 


1.52 


0.55 


1.65 


0.38 


0.13 


0.07 


15.09 


0.00 


0.00 


0.91 


CJS3-5 


0.42 


13.39 


33.61 


6.95 


25.59 


1.03 


1.49 


0.55 


1.62 


0.37 


0.13 


0.07 


14.79 


0.00 


0.00 


0.91 


CJS3-6 


1.27 


12.10 


38.97 


6.46 


23.11 


0.93 


1.35 


0.50 


1.46 


0.33 


0.11 


0.07 


13.36 


0.00 


0.00 


0.91 


CJS4-4 


0.08 


14.28 


31.63 


5.36 


27.30 


1.09 


1.59 


0.58 


1.73 


0.39 


0.13 


0.08 


15.77 


0.00 


0.00 


0.91 


CJS4.5 


0.55 


13.05 


34.98 


6.96 


24.93 


1.00 


1.46 


0.54 


1.58 


0.36 


0.12 


0.07 


14.41 


0.00 


0.00 


0.91 


EJSl-4 


0.00 


14.33 


30.48 


6.39 


27.38 


1.10 


1.60 


0..58 


1.73 


0.40 


0.13 


0.08 


15.82 


0.00 


0.00 


0.91 


EJSl-5 


0.07 


14.24 


31.90 


5.29 


27.20 


1.09 


1.59 


0.58 


1.72 


0.39 


0.13 


0.08 


15.73 


0.00 


0.00 


0.91 


EJSl-6 


0.01 


13.74 


32.09 


7.35 


26.26 


1.05 


1.53 


0.57 


1.66 


0.38 


0.13 


0.07 


15.17 


0.00 


0.00 


0.91 


EJS2-4 


0.04 


14.55 


31.86 


4.00 


27.81 


1.11 


1.62 


0.58 


1.76 


0.40 


0.13 


0.08 


16.07 


0.00 


0.00 


0.91 


EJS2-5 


0.00 


14.39 


30.94 


5.69 


27.47 


1.10 


1.60 


0.59 


1.74 


0.40 


0.13 


0.08 


15.89 


0.00 


0.00 


0.91 


EJS2-6 


0.00 


14.00 


30.97 


7.30 


26.77 


1.07 


1.56 


0.58 


1.70 


0.39 


0.13 


0.07 


15.47 


0.00 


0.00 


0.91 
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Table 8 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




EJS2-7 


0.67 


12.81 


36.09 


6.86 


24.42 


0.98 


1.43 


0.53 


1.55 


0.35 


0.12 


0.07 


14.14 


0.00 


0.00 


0.91 


EJS3-4 


0.00 


14.46 


31.22 


5.04 


27.65 


1.11 


1.61 


0.58 


1.75 


0.40 


0.13 


0.08 


15.98 


0.00 


0.00 


0.91 


EJS3-5 


0.00 


14.14 


30.41 


7.25 


27.04 


1.08 


1.58 


0.58 


1.71 


0.39 


0.13 


0.08 


15.62 


0.00 


0.00 


0.91 


EJS3-6 


0.00 


14.60 


32.95 


1.81 


25.26 


2.41 


3.44 


0.50 


1.61 


0.36 


0.12 


0.16 


16.81 


0.00 


0.00 


0.87 


EJS3-7 


4.74 


7.89 


56.34 


4.12 


15.10 


0.61 


0.88 


0.33 


0.96 


0.22 


0.08 


0.05 


8.72 


0.00 


0.00 


0.91 


EJS4-4 


0.00 


14.69 


30.95 


4.32 


28.07 


1.13 


1.64 


0.59 


1.78 


0.41 


0.14 


0.08 


16.23 


0.00 


0.00 


0.91 


EJS4-5 


0.00 


14.36 


30.33 


6.37 


27.46 


1.10 


1.60 


0..58 


1.74 


0.40 


0.13 


0.08 


15.87 


0.00 


0.00 


0.91 


EJS4-6 


0.02 


13.73 


32.10 


7.33 


26.26 


1.05 


1.53 


0.57 


1.66 


0.38 


0.13 


0.07 


15.17 


0.00 


0.00 


0.91 


EJS4.7 


0.01 


13.82 


31.81 


7.30 


26.40 


1.06 


1.54 


0.57 


1.67 


0.38 


0.13 


0.07 


15.26 


0.00 


0.00 


0.91 


t=1.5xl0'' years 


CJSl-4 


1.23 


12.56 


37.60 


5.84 


23.99 


0.96 


1.40 


0.52 


1.52 


0.35 


0.14 


0.06 


13.87 


0.00 


0.00 


0.91 


CJSl-5 


0.94 


12.64 


36.65 


6.73 


24.15 


0.97 


1.41 


0.52 


1.53 


0.35 


0.11 


0.06 


13.96 


0.00 


0.00 


0.91 


CJSl-6 


1.58 


11.72 


40.49 


6.29 


22.39 


0.90 


1.31 


0.48 


1.42 


0.32 


0.13 


0.06 


12.94 


0.00 


0.00 


0.91 


CJS2-4 


0.15 


14.23 


31.61 


5.52 


27.20 


1.09 


1.59 


0.59 


1.72 


0.39 


0.12 


0.07 


15.72 


0.00 


0.00 


0.91 


CJS2-5 


0.26 


13.45 


33.32 


7.18 


25.68 


1.03 


1.50 


0.55 


1.63 


0.37 


0.12 


0.06 


14.85 


0.00 


0.00 


0.91 


CJS2-6 


1.15 


12.41 


37.59 


6.57 


23.72 


0.95 


1.38 


0.51 


1.50 


0.34 


0.11 


0.06 


13.71 


0.00 


0.00 


0.91 


CJS2-7 


1.72 


11.52 


41.35 


6.18 


22.00 


0.88 


1.28 


0.48 


1.39 


0.32 


0.10 


0.05 


12.72 


0.00 


0.00 


0.91 


CJS3-4 


0.97 


13.03 


35.94 


5.67 


24.90 


1.00 


1.45 


0.54 


1.58 


0.36 


0.11 


0.06 


14.39 


0.00 


0.00 


0.91 


CJS3-5 


0.67 


12.91 


35.53 


6.91 


24.67 


0.99 


1.44 


0.53 


1.56 


0.36 


0.11 


0.06 


14.26 


0.00 


0.00 


0.91 


CJS3-6 


3.16 


9.76 


48.57 


5.24 


18.66 


0.75 


1.09 


0.40 


1.18 


0.27 


0.08 


0.05 


10.78 


0.00 


0.00 


0.91 


CJS4-4 


0.25 


13.74 


32.26 


6.95 


26.25 


1.05 


1.53 


0.57 


1.66 


0.38 


0.12 


0.06 


15.17 


0.00 


0.00 


0.91 


CJS4-5 


1.50 


11.82 


40.07 


6.33 


22.59 


0.91 


1.32 


0.49 


1.43 


0.33 


0.10 


0.05 


13.06 


0.00 


0.00 


0.91 


EJSl-4 


0.04 


13.92 


31.31 


7.28 


26.62 


1.07 


1.55 


0.57 


1.69 


0.38 


0.14 


0.07 


15.38 


0.00 


0.00 


0.91 


EJSl-5 


0.60 


13.29 


34.63 


6.22 


25.39 


1.02 


1.48 


0.55 


1.61 


0.37 


0.14 


0.07 


14.67 


0.00 


0.00 


0.91 


EJSl-6 


0.33 


13.22 


34.29 


7.10 


25.27 


1.01 


1.48 


0.55 


1.60 


0.37 


0.14 


0.06 


14.60 


0.00 


0.00 


0.91 
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Table 8 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




EJS2-4 


0.05 


14.22 


31.92 


5.38 


27.16 


1.09 


1.59 


0.58 


1.72 


0.39 


0.14 


0.07 


15.71 


0.00 


0.00 


0.91 


EJS2-5 


0.03 


13.95 


31.33 


7.16 


26.67 


1.07 


1.56 


0.58 


1.69 


0.39 


0.14 


0.07 


15.41 


0.00 


0.00 


0.91 


E.JS2-6 


0.04 


13.73 


32.14 


7.34 


26.23 


1.05 


1.53 


0.57 


1.66 


0..38 


0.14 


0.06 


15.16 


0.00 


0.00 


0.91 




O.OZ 


"7 1 Q 

/.io 


oy.oy 


o.o4 


1 Q fit; 
io.DO 


U.OO 


U.oU 


0.30 


0.87 


n on 


U. io 


U.UD 


I .OO 


U.UU 


n nn 
U.UU 


n nn 
u.yU 


EJS3-4 


0.05 


14.22 


31.11 


6.20 


27.16 


1.09 


1.59 


0.59 


1.72 


0.39 


0.14 


0.07 


15.70 


0.00 


0.00 


0.91 


EJS3-5 


0.04 


13.72 


32.16 


7.32 


26.23 


1.05 


1.53 


0.57 


1.66 


0.38 


0.14 


0.07 


15.16 


0.00 


0.00 


0.91 


EJS3-6 


0.59 


13.04 


35.04 


6.92 


24.91 


1.00 


1.45 


0.54 


1.58 


0.36 


0.13 


0.15 


14.40 


0.00 


0.00 


0.91 




o.uy 


'7/1/1 
/ .44 


oo. Id 


Q OQ 

o.yo 


1 A on 
14. zU 


U.o ( 


U.OO 


0.31 


0.90 


n oi 


n nn 
U.Uy 


n c\A 
U.U4 


Q OO 


n nn 
U.UU 


n nn 
U.UU 


n ni 
U.yl 


EJS4-4 


0.02 


14.22 


30.97 


6.36 


27.16 


1.09 


1.59 


0.59 


1.72 


0.39 


0.15 


0.07 


15.70 


0.00 


0.00 


0.91 


EJS4-5 


0.01 


14.02 


31.05 


7.14 


26.81 


1.07 


1.56 


0.58 


1.70 


0.39 


0.12 


0.06 


15.49 


0.00 


0.00 


0.91 




n 1 n 
U. 19 




OO. lo 




Zd. i 4 


i.Uo 


1 

i . DU 


0.56 


1.63 


U.o i 


n 1 /I 
U. 14 


U.UD 


1 A Qn 
14. oy 


n nn 
U.UU 


n nn 
U.UU 


n ni 
U.yl 


E.JSl-7 


0.27 


13.35 


33.71 


7.16 


25.52 


1.02 


1.49 


0.55 


1.62 


0.37 


0.14 


0.06 


14.75 


0.00 


0.00 


0.91 
















t.=2x 


10"' yc.-ars 
















CJSl-4 


1.65 


11.86 


39.98 


6.11 


22.66 


0.91 


1.32 


0.49 


1.44 


0.33 


0.11 


0.04 


13.10 


0.00 


0.00 


0.91 


CJSl-5 


1.09 


12.39 


37.74 


6.62 


23.65 


0.95 


1.38 


0.51 


1.50 


0.34 


0.11 


0.06 


13.68 


0.00 


0.00 


0.91 


CJSl-6 


2.03 


11.15 


42.88 


5.98 


21.29 


0.85 


1.24 


0.46 


1.35 


0.31 


0.10 


0.05 


12.31 


0.00 


0.00 


0.91 


CJS2-4 


0.35 


13.70 


32.53 


6.77 


26.17 


1.05 


1.53 


0.57 


1.66 


0.38 


0.13 


0.05 


15.13 


0.00 


0.00 


0.91 


CJS2-5 


0.62 


12.92 


35.49 


6.94 


24.70 


0.99 


1.44 


0.53 


1.56 


0.36 


0.12 


0.05 


14.27 


0.00 


0.00 


0.91 


CJS2-6 


1.30 


12.13 


38.79 


6.49 


23.16 


0.93 


1.35 


0.50 


1.47 


0.33 


0.10 


0.06 


13.39 


0.00 


0.00 


0.91 


CJS2-7 


3.16 


9.78 


48.51 


5.25 


18.67 


0.75 


1.09 


0.40 


1.18 


0.27 


0.08 


0.05 


10.80 


0.00 


0.00 


0.91 


CJS3-4 


1.16 


12.56 


37.13 


6.38 


23.99 


0.96 


1.40 


0.52 


1.52 


0.35 


0.12 


0.05 


13.87 


0.00 


0.00 


0.91 


CJS3-5 


1.20 


12.21 


38.49 


6.55 


23.31 


0.93 


1.36 


0.50 


1.48 


0.34 


0.11 


0.06 


13.48 


0.00 


0.00 


0.91 


CJS3-6 


4.35 


8.36 


54.35 


4.48 


15.95 


0.64 


0.93 


0.34 


1.01 


0.23 


0.07 


0.04 


9.23 


0.00 


0.00 


0.91 


CJS4-4 


0.77 


12.91 


35.63 


6.72 


24.66 


0.99 


1.44 


0.53 


1.56 


0.36 


0.11 


0.06 


14.26 


0.00 


0.00 


0.91 


CJS4.5 


2.47 


10.65 


44.92 


5.71 


20.33 


0.82 


1.19 


0.44 


1.29 


0.29 


0.09 


0.05 


11.76 


0.00 


0.00 


0.91 


EJSl-4 


0.13 


13.62 


32.67 


7.22 


26.01 


1.04 


1.52 


0.56 


1.65 


0.38 


0.13 


0.05 


15.04 


0.00 


0.00 


0.91 
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Table 8 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




EJSl-5 


0.75 


12.94 


35.50 


6.73 


24.72 


0.99 


1.44 


0.53 


1.57 


0.36 


0.12 


0.05 


14.29 


0.00 


0.00 


0.91 


EJSl-6 


2.40 


10.73 


44.61 


5.75 


20.48 


0.82 


1.20 


0.44 


1.30 


0.30 


0.09 


0.05 


11.84 


0.00 


0.00 


0.91 


EJS2-4 


0.27 


13.69 


32.75 


6.67 


26.15 


1.05 


1.53 


0.57 


1.66 


0.38 


0.13 


0.05 


15.12 


0.00 


0.00 


0.91 


EJS2-5 


0.19 


13.60 


32.70 


7.22 


25.97 


1.04 


1.52 


0.56 


1.65 


0.38 


0.13 


0.05 


15.01 


0.00 


0.00 


0.91 


EJS2-6 


0.24 


13.39 


33.62 


7.18 


25.56 


1.02 


1.49 


0.55 


1.62 


0.37 


0.12 


0.06 


14.78 


0.00 


0.00 


0.91 


EJS2-7 


5.57 


6.84 


60.65 


3.67 


13.05 


0.52 


0.76 


0.28 


0.83 


0.19 


0.06 


0.03 


7.55 


0.00 


0.00 


0.91 


EJS3-4 


0.13 


13.80 


32.06 


7.00 


26.37 


1.06 


1.54 


0.57 


1.67 


0.38 


0.13 


0.06 


15.24 


0.00 


0.00 


0.91 


EJS3-5 


0.13 


13.57 


32.86 


7.27 


25.90 


1.04 


1.51 


0.56 


1.64 


0.37 


0.13 


0.06 


14.98 


0.00 


0.00 


0.91 


EJS3-6 


0.20 


13.43 


33.39 


7.21 


25.67 


1.03 


1.50 


0.55 


1.63 


0.37 


0.12 


0.06 


14.84 


0.00 


0.00 


0.91 


E.JS3-7 


0.06 


8.37 


30.35 


1.49 


32.96 


3.49 


2.55 


0.97 


0.01 


0.30 


0.12 


4.53 


13.99 


0.00 


0.79 


0.60 


EJS4-4 


0.07 


13.90 


31.57 


7.12 


26.56 


1.06 


1.55 


0.57 


1.68 


0.38 


0.13 


0.05 


15.35 


0.00 


0.00 


0.91 


EJS4-5 


0.04 


13.74 


32.24 


7.20 


26.24 


1.05 


1.53 


0.57 


1.66 


0..38 


0.13 


0.05 


15.17 


0.00 


0.00 


0.91 


EJSl-6 


0.25 


13.29 


31.10 


7.13 


25.36 


1.02 


1.18 


0.55 


1.61 


0.37 


0.11 


0.06 


11.67 


0.00 


0.00 


0.91 


EJSi-7 


l.:-!7 


11.98 


39. 10 


(). 13 


22.89 


0.92 


1.:-!! 


0. 19 


1.15 


0.3:i 


0.10 


0.0() 


i:-!.2:i 


0.00 


0.00 


0.91 


t=2.5xl0" years 


CJSl-4 


1.88 


11.52 


41.25 


6.11 


22.01 


0.88 


1.29 


0.48 


1.39 


0.32 


0.10 


0.05 


12.72 


0.00 


0.00 


0.91 


CJSl-5 


1.31 


12.06 


39.08 


6.46 


23.04 


0.92 


1.35 


0.50 


1.46 


0.33 


0.10 


0.06 


13.32 


0.00 


0.00 


0.91 


CJSl-6 


2.82 


10.21 


46.72 


5.48 


19.50 


0.78 


1.14 


0.42 


1.24 


0.28 


0.09 


0.05 


11.28 


0.00 


0.00 


0.91 


CJS2-4 


0.47 


13.35 


33.67 


7.03 


25.52 


1.02 


1.49 


0.55 


1.62 


0.37 


0.12 


0.06 


14.75 


0.00 


0.00 


0.91 


CJS2-5 


0.88 


12.60 


36.88 


6.76 


24.05 


0.96 


1.40 


0.52 


1.52 


0.35 


0.11 


0.06 


13.91 


0.00 


0.00 


0.91 


CJS2-6 


1.86 


11.41 


41.73 


6.12 


21.80 


0.87 


1.27 


0.47 


1.38 


0.32 


0.10 


0.05 


12.60 


0.00 


0.00 


0.91 


CJS2-7 


4.01 


8.75 


52.73 


4.70 


16.71 


0.67 


0.98 


0.36 


1.06 


0.24 


0.08 


0.04 


9.67 


0.00 


0.00 


0.91 


CJS3-4 


1.24 


12.34 


37.83 


6.54 


23.59 


0.95 


1.38 


0.51 


1.49 


0..34 


0.11 


0.06 


13.63 


0.00 


0.00 


0.91 


CJS3-5 


2.00 


11.22 


42.54 


6.02 


21.43 


0.86 


1.25 


0.46 


1.36 


0.31 


0.10 


0.05 


12.39 


0.00 


0.00 


0.91 


CJS3-6 


4.47 


8.21 


54.96 


4.40 


15.68 


0.63 


0.92 


0.34 


0.99 


0.23 


0.07 


0.04 


9.07 


0.00 


0.00 


0.91 


CJS4-4 


0.93 


12.63 


36.69 


6.73 


24.13 


0.97 


1.41 


0.52 


1.53 


0.35 


0.11 


0.06 


13.95 


0.00 


0.00 


0.91 
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Table 8 - continued from previous page 



Simulation 


H 


Mg 


O 


S 


Fe 


Al 


Ca 


Na 


Ni 


Cr 


P 


Ti 


Si 


N 


c 


Mg/Si 




wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 


wt% 




CJS4-5 


3.40 


9.49 


49.71 


5.09 


18.12 


0.73 


1.06 


0.39 


1.15 


0.26 


0.08 


0.04 


10.48 


0.00 


0.00 


0.91 


EJSl-4 


0.30 


13.35 


33.74 


7.13 


25.50 


1.02 


1.49 


0.55 


1.62 


0.37 


0.12 


0.06 


14.74 


0.00 


0.00 


0.91 


EJSl-5 


0.96 


12.58 


36.87 


6.72 


24.04 


0.96 


1.40 


0.52 


1.52 


0.35 


0.11 


0.06 


13.90 


0.00 


0.00 


0.91 


EJSl-6 


3.51 


9.34 


50.34 


5.01 


17.83 


0.72 


1.04 


0.39 


1.13 


0.26 


0.08 


0.04 


10.32 


0.00 


0.00 


0.91 


E.JS2-4 


0.36 


13.46 


33.28 


7.04 


25.73 


1.03 


1.50 


0.56 


1.63 


0.37 


0.12 


0.06 


14.86 


0.00 


0.00 


0.91 


EJS2-5 


0.80 


12.74 


36.25 


6.82 


24.34 


0.98 


1.42 


0.53 


1.54 


0.35 


0.11 


0.06 


14.07 


0.00 


0.00 


0.91 


EJS2-6 


0.52 


13.01 


35.20 


6.98 


24.84 


1.00 


1.45 


0.54 


1.57 


0.36 


0.11 


0.06 


14.37 


0.00 


0.00 


0.91 


EJS2-7 


5.57 


6.84 


60.65 


3.67 


13.05 


0.52 


0.76 


0.28 


0.83 


0.19 


0.06 


0.03 


7.55 


0.00 


0.00 


0.91 


E.JS3-4 


0.22 


13.58 


32.71 


7.20 


25.96 


1.04 


1.52 


0.56 


1.64 


0.38 


0.12 


0.06 


15.00 


0.00 


0.00 


0.91 


EJS3-5 


0.21 


13.42 


33.50 


7.19 


25.61 


1.03 


1.50 


0.55 


1.62 


0.37 


0.12 


0.06 


14.81 


0.00 


0.00 


0.91 


EJS3-6 


1.66 


11.61 


40.94 


6.23 


22.19 


0.89 


1.30 


0.48 


1.41 


0.32 


0.10 


0.05 


12.82 


0.00 


0.00 


0.91 


EJS3-7 


5.35 


7.29 


58.65 


3.91 


13.91 


0.56 


0.81 


0.30 


0.88 


0.20 


0.06 


0.03 


8.05 


0.00 


0.00 


0.91 


EJS4^4 


0.15 


13.65 


32.44 


7.25 


26.09 


1.05 


1.52 


0.56 


1.65 


0.38 


0.12 


0.06 


15.08 


0.00 


0.00 


0.91 


EJS4.5 


0.12 


13.60 


32.75 


7.24 


25.96 


1.04 


1.52 


0.56 


1.65 


0.38 


0.12 


0.06 


15.01 


0.00 


0.00 


0.91 


EJS4-6 


0.83 


12.61 


36.88 


6.76 


24.06 


0.96 


1.41 


0.52 


1.52 


0.35 


0.11 


0.06 


13.92 


0.00 


0.00 


0.91 


EJS4-7 


2.54 


10.53 


45.40 


5.65 


20.12 


0.81 


1.17 


0.43 


1.27 


0.29 


0.09 


0.05 


11.63 


0.00 


0.00 


0.91 


t=3xlO'' years 


CJSl-4 


1.99 


11.36 


41.91 


6.07 


21.69 


0.87 


1.27 


0.47 


1.37 


0.31 


0.10 


0.05 


12.54 


0.00 


0.00 


0.91 


CJSl-5 


1.70 


11.60 


41.00 


6.22 


22.14 


0.89 


1.29 


0.48 


1.40 


0.32 


0.10 


0.05 


12.81 


0.00 


0.00 


0.91 


CJSl-6 


4.20 


8.53 


53.61 


4.58 


16.31 


0.65 


0.95 


0.35 


1.03 


0.24 


0.07 


0.04 


9.43 


0.00 


0.00 


0.91 


CJS2-4 


0.57 


13.15 


34.43 


7.03 


25.14 


1.01 


1.47 


0.54 


1.59 


0.36 


0.11 


0.06 


14.53 


0.00 


0.00 


0.91 


CJS2-5 


2.11 


11.10 


43.08 


5.95 


21.18 


0.85 


1.24 


0.46 


1.34 


0.31 


0.10 


0.05 


12.25 


0.00 


0.00 


0.91 


CJS2-6 


2.04 


11.18 


42.71 


6.00 


21.35 


0.86 


1.25 


0.46 


1.35 


0.31 


0.10 


0.05 


12.35 


0.00 


0.00 


0.91 


CJS2-7 


4.86 


7.76 


56.78 


4.17 


14.83 


0.59 


0.87 


0.32 


0.94 


0.21 


0.07 


0.04 


8.57 


0.00 


0.00 


0.91 


CJS3-4 


0.62 


12.96 


35.34 


6.93 


24.75 


0.99 


1.45 


0.53 


1.57 


0.36 


0.11 


0.06 


14.31 


0.00 


0.00 


0.91 


CJS3-5 


2.55 


10.55 


45.31 


5.66 


20.15 


0.81 


1.18 


0.44 


1.28 


0.29 


0.09 


0.05 


11.65 


0.00 


0.00 


0.91 
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Simulation 


1 1 

wt% 


Mg 
wt% 


° t 


wt% 


re 
wt% 


A 1 

Wt% 


L^a 
wt% 


IN a 
wt% 


iNl 

Wt% 


wt% 


T 1 

ir 
wt% 


Ti 

wt% 


CI 

»1 

wt% 


TVT 
IN 

Wt% 


e 1 

wt% 


Mg/bi 


CJS3-6 


5.47 


7.05 


59.67 


3.78 


13.47 


0.54 


0.79 


0.29 


0.85 


0.19 


0.06 


0.03 


7.79 


0.00 


0.00 


0.91 


CJS4-4 


1.13 


12.35 


37.86 


6.62 


23.58 


0.95 


1.38 


0.51 


1.49 


0..34 


0.11 


0.06 


13.63 


0.00 


0.00 


0.91 


CJS4-5 


4.26 


8.47 


53.88 


4.54 


16.17 


0.65 


0.94 


0.35 


1.02 


0.23 


0.07 


0.04 


9.35 


0.00 


0.00 


0.91 


EJSl-4 


0.49 


13.08 


34.91 


7.00 


24.96 


1.00 


1.46 


0.54 


1.58 


0.36 


0.11 


0.06 


14.44 


0.00 


0.00 


0.91 


EJSl-5 


1.66 


11.70 


40.49 


6.28 


22.36 


0.90 


1.31 


0.48 


1.42 


0.32 


0.10 


0.05 


12.93 


0.00 


0.00 


0.91 


EJSl-6 


2.01 


11.15 


42.85 


5.99 


21.31 


0.85 


1.24 


0.46 


1.35 


0.31 


0.10 


0.05 


12.32 


0.00 


0.00 


0.91 


EJS2-4 


0.31 


13.38 


33.61 


7.13 


25.56 


1.02 


1.49 


0.55 


1.62 


0.37 


0.12 


0.06 


14.78 


0.00 


0.00 


0.91 


EJS2-5 


0.46 


13.19 


34.38 


7.06 


25.18 


1.01 


1.47 


0.54 


1.60 


0.36 


0.11 


0.06 


14.57 


0.00 


0.00 


0.91 


EJS2-6 


1.31 


12.03 


39.25 


6.45 


22.97 


0.92 


1.34 


0.50 


1.46 


0.33 


0.10 


0.06 


13.28 


0.00 


0.00 


0.91 


EJS2-7 


5.57 


6.84 


60.63 


3.67 


13.06 


0.52 


0.76 


0.28 


0.83 


0.19 


0.06 


0.03 


7.55 


0.00 


0.00 


0.91 


EJS3-4 


0.44 


13.27 


34.04 


7.10 


25.33 


1.02 


1.48 


0.55 


1.60 


0.37 


0.11 


0.06 


14.65 


0.00 


0.00 


0.91 


EJS3-5 


0.50 


13.03 


35.11 


6.99 


24.88 


1.00 


1.45 


0.54 


1.58 


0.36 


0.11 


0.06 


14.39 


0.00 


0.00 


0.91 


EJS3-6 


2.88 


10.12 


47.08 


5.44 


19.34 


0.78 


1.13 


0.42 


1.23 


0.28 


0.09 


0.05 


11.18 


0.00 


0.00 


0.91 


EJS3-7 


5.37 


7.27 


58.74 


3.90 


13.86 


0.56 


0.81 


0.30 


0.88 


0.20 


0.06 


0.03 


8.02 


0.00 


0.00 


0.91 


EJS4^4 


0.25 


13.47 


33.18 


7.21 


25.73 


1.03 


1.50 


0.56 


1.63 


0.37 


0.12 


0.06 


14.88 


0.00 


0.00 


0.91 


EJS4-5 


0.36 


13.28 


34.06 


7.10 


25.34 


1.02 


1.48 


0.55 


1.61 


0.37 


0.11 


0.06 


14.66 


0.00 


0.00 


0.91 


EJS4-6 


1.23 


12.09 


38.97 


6.50 


23.13 


0.93 


1.35 


0.50 


1.47 


0.33 


0.10 


0.06 


13.35 


0.00 


0.00 


0.91 


EJS4-7 


3.20 


9.72 


48.78 


5.21 


18.55 


0.74 


1.08 


0.40 


1.18 


0.27 


0.08 


0.04 


10.73 


0.00 


0.00 


0.91 
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Figure 8: Normalized planetary abundances for CJS2 simulated terrestrial planets. Values 
are shown for each of seven time steps considered with the following color scheme: black = 
2.5x10^ years, red = 5x10^ years, green = 1x10® years, pink = 1.5x10® years, light blue = 
2x10® years, yellow = 2.5x10® years and dark blue = 3x10® years. 
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Figure 9: Normalized planetary abundances for CJS3 and CJS4 simulated terrestrial planets. 
Values are shown for each of seven time steps considered with the following color scheme: 
black = 2.5x10^ years, red = 5x10^ years, green = 1x10^ years, pink = 1.5x10^ years, 
hght blue = 2x10^ years, yellow = 2.5x10^ years and dark blue = 3x10^ years. Left: CJS3 
simulation results. Right: CJS4 simulation results. 
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Figure 10: Normalized planetary abundances for EJS2 simulated terrestrial planets. Values 
are shown for each of seven time steps considered with the following color scheme: black = 
2.5x10^ years, red = 5x10^ years, green = 1x10^ years, pink = 1.5x10® years, light blue = 
2x10® years, yellow = 2.5x10® years and dark blue = 3x10® years. 
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Figure 11: Normalized planetary abundances for EJS3 simulated terrestrial planets. Values 
are shown for each of seven time steps considered with the following color scheme: black = 
2.5x10^ years, red = 5x10^ years, green = 1x10® years, pink = 1.5x10® years, light blue = 
2x10® years, yellow = 2.5x10® years and dark blue = 3x10® years. 
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Figure 12: Normalized planetary abundances for EJS4 simulated terrestrial planets. Values 
are shown for each of seven time steps considered with the following color scheme: black = 
2.5x10^ years, red = 5x10^ years, green = 1x10® years, pink = 1.5x10'' years, light blue = 
2x10® years, yellow = 2.5x10® years and dark blue = 3x10® years. 
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